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Abstract 

In this second paper we define a Post-Minkowskian (PM) weak field approximation leading to 
a linearization of the Hamilton equations of ADM tetrad gravity in the York canonical basis in a 
family of non-harmonic 3-orthogonal Schwinger time gauges. The York time 3 K (the relativistic 
inertial gauge variable, not existing in Newtonian gravity, parametrizing the family and connected 
to the freedom in clock synchronization, i.e. to the definition of the instantaneous 3-spaces) is 
put equal to an arbitrary numerical function. The matter are point particles, with a Grassmann 
regularization of self-energies, and the electro-magnetic field in the radiation gauge: a ultraviolet 
cutoff allows a consistent linearization, which is shown to be the lowest order of a Hamiltonian 
Post-Minkowskian (HPM) expansion. 

We solve the constraints and the Hamilton equations for the tidal variables and we find Post- 
Minkowskian gravitational waves with asymptotic background (and the correct quadrupole emis- 
sion formula) propagating on dynamically determined non-Euclidean 3-spaces. The conserved 
ADM energy and the Grassmann regularizzation of self-energies imply the correct energy balance. 
A generalized transverse-traceless gauge can be identified and the main tools for the detection 
of gravitational waves are reproduced in these non-harmonic gauges. In conclusion we get a PM 
solution for the gravitational field and we identify a class of PM Einstein space-times, which will be 
studied in more detail in a third paper together with the PM equations of motion for the particles 
and their Post-Newtonian expansion (but in absence of the electro-magnetic field) . 

Finally we make a discussion on the gauge problem in general relativity to understand which 
type of experimental observations may lead to a preferred choice for the inertial gauge variable 
S K in the PM space-times. In the third paper we will show that this choice is connected with the 
problem of dark matter. 
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I. INTRODUCTION 



In this paper we will study the linearization of the Hamiltonian formulation of ADM 
tetrad gravity given in Ref . [1] (quoted as paper I in what follows) to get a Post-Minkowskian 
(PM) description of gravitational waves (GW) in non-harmonic gauges in asymptotically 
Minkowskian space-times with the asymptotic Minkowski metric used as a background. We 
are able to reproduce all the main properties of GW's, which are usually derived in the 
standard Lagrangian approach to GW's in harmonic gauges with a Post-Newtonian (PN) 
expansion (see Appendix A for a review). 

We define a new Hamiltonian Post-Minkowskiam (HPM) approach to the description 
of GW's in which we do not assume the decomposition = 4 rj^ u + 4 /i At „, but we use 
an asymptotic Minkowski background 4-metric A rj^ asyrn ) at spatial infinity in a certain 
family of asymptotically Minkowskian space-times. As shown in paper I, this approach is 
based on ADM tetrad gravity in these space-times followed by a canonical transformation 
to a York canonical basis, which diagonalizes the York-Lichnerowitz formulation of general 
relativity [2, 3] allowing a clean separation between physical tidal degrees of freedom and 
inertial gauge variables inside the gravitational field in a 3-covariant way inside the non- 
Euclidean instantaneous 3-spaces (the Cauchy surfaces for the tidal variables): while in 
special relativity its shape as a sub-manifold of Minkowski space-time is arbitrary [4] (the 
gauge freedom in the choice of the convention for clock synchronization), in general relativity 
the 3-space is dynamically determined [5] except for the trace of its extrinsic curvature, the 
inertial gauge variable called York time (it describes the general relativistic remnant of the 
gauge freedom in clock synchronization). 

In paper I we have explicitly evaluated the Hamilton equations of ADM tetrad gravity 
coupled to dynamical matter in the Schwinger time gauges in the York canonical basis defined 
in Ref. [6]. This paper was the final development of previous researches [7-10] in canonical 
gravity. The matter consists of the electro-magnetic field in the radiation gauge and of N 
dynamical (not test) massive point particles 1 . The Grassmann- valued electric charges and 
signs of the energy regularize both the electro-magnetic and gravitational self-energies (see 
Refs.[4, 12-15] for the electro-magnetic case in special relativity) in the equations of motion 
for the N-body problem avoiding both the gravitational and electro-magnetic self-energies. 
As a consequence the action principle for the particles given in paper I is well posed (there 
are no essential singularities on the particle world-lines). 

Dirac theory of constraints is taken into account at every step, in particular when making 
gauge-fixings (instead in the standard ADM approach and numerical gravity the gauge fix- 
ings are chosen only on the basis of convenience). This will allow to make a clear distinction 
between the instantaneous action-at-a-distance effects implied by the super-Hamiltonian 
and super-momentum constraints (like it happens in every gauge theory: with the electro- 
magnetic field in the radiation gauge this is implied by the Gauss law constraint in presence 



In a future paper we will try to describe compact extended objects and their self-gravity (see footnote 16 
of Appendix A) starting from balls of perfect fluids, whose Lagrangian description will be the extension 
to tetrad gravity of the special relativistic one given in Ref[ll]. The Grassmann regularization is a semi- 
classical way out from causality problems like the ones arising when considering the point limit of extended 
models for the electron in classical electro-dynamics. 
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of matter) and retarded tidal effects. 

We use a family of non-harmonic 3-orthogonal Schwinger time-gauges in which the gauge 
fixings imply only elliptic-type equations at a given time for the gauge variables. As a con- 
sequence, we avoid the wave equations for gauge variables like the lapse and shift functions 
present in the harmonic gauges and requiring initial data in the asymptotic past. 

ADM tetrad gravity is formulated in an arbitrary admissible 3+1 splitting of the globally 
hyperbolic space-time, i.e. in a foliation with instantaneous space-like 3-spaces tending to 
a Minkowski space-like hyper-plane at spatial infinity in a direction independent way: they 
correspond to a clock synchronization convention and each one of them can be used as a 
Cauchy surface for field equations. As shown in Refs.[l, 7] the absence of super-translations 
implies that the SPI group of asymptotic symmetries is reduced to the asymptotic ADM 
Poincare' group 2 and the allowed 3+1 splittings must have the instantaneous 3-spaces 
tending to asymptotic special-relativistic Wigner hyper-planes orthogonal to the ADM 4- 
momentum in a direction-independent way (see Ref. [4] for these non-inertial rest frames in 
special relativity). At spatial infinity there are asymptotic inertial observers, carrying a flat 
tetrad { A V^ e A e B = 4 Vab, e M Q( g 7 e" = e£), whose spatial axes can be identified with 
the fixed stars of star catalogues. 

We use radar 4-coordinates a A = (a T = r; a r ), A = r,r, adapted to the admissible 3+1 
splitting of the space-time and centered on an arbitrary time-like observer a+(r) (origin of 
the 3-coordinates a r ): they define a non-inertial frame centered on the observer, so that 
they are observer and frame- dependent. The time variable r is an arbitrary monotonically 
increasing function of the proper time given by the atomic clock carried by the observer. The 
instantaneous 3-spaces identified by this convention for clock synchronization are denoted 
S T . The transformation u A i— > a+ = z^{r, a r ) to world 4-coordinates defines the embedding 
z^(t, c?) of the Riemannian instantaneous 3-spaces S T into the space-time. By choosing world 
4-coordinates centered on the time-like observer, whose world-line is the time axis, we have 
a+(r) = (x°(t);0): the condition x°(t) = const, is equivalent to t — const, and identifies 
the instantaneous 3-space S T . If the time-like observer coincides with an asymptotic inertial 
observer x^(r) = x^ + e^r with = (1;0), e£f = (0;8 l r ), x% = (x°;0), then the natural 
embedding describing the given 3+1 splitting of space-time is z^{t, a r ) = x^ + a A and the 
world 4- metric is = e A a qab (e^ are flat asymptotic cotetrads, e A e^ B = 8 A , e A e y A = 8%). 

From now on we shall denote the curvilinear 3-coordinates a r with the notation a for the 
sake of simplicity. Usually the convention of sum over repeated indices is used, except when 
there are too many summations. 

The 4-metric a Qab has signature e (H ) with e = ± (the particle physics, e = +, and 

general relativity, e = — , conventions). Flat indices (a), a = o,a, are raised and lowered by 
the flat Minkowski metric 4 ^( Q )(/3) = e (H ). We define 4 ^( a )(6) — — e ^(a)(6) with a positive- 
definite Euclidean 3-metric. On each instantaneous 3-space S T we have that the 4-metric has 
a direction-independent limit to the flat Minkowski 4-metric (the asymptotic background) 
at spatial infinity A gAB(j,a) — > 4 i]AB(asym) — e(H )• This asymptotic 4-metric allows 



2 It reduces to the special relativistic Poincare' group of the given matter in non-inertial frames of Minkowski 
space-time when the Newton constant G is switched off. 
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to define both a flat d'Alambertian □ = — A and a flat Laplacian A = d% on S r 
{pA = q~k)- We will also need the flat distribution c(a,a) = ^ 5 3 (<x, a') = — \g-ff'\ with 
\a — a\ = (°"" ~~ cr ' u ) 2 i where 5 3 (<j, a) is the Dirac delta function on the 3-manifold 

In this paper we study the linearization of the Hamilton equations of ADM tetrad gravity 
in the York canonical basis with the weak field approximation in the family of non-harmonic 
3-orthogonal gauges in which the gauge variable (a relativistic inertial effect) York time 
3 K(t, a) is an arbitrary numerical function and the 3-coordinates in £ r are chosen so that 
the 3-metric is everywhere diagonal. This will be a starting point before facing higher orders 
in the HPM expansion defined at the end of Section III (at higher orders regular izat ion 
problems may arise). 

Since the gravitational gauge freedom is fixed, we get well defined Hamilton equations 
for the matter. In particular this avoids the introduction of "ad hoc" Lagrangians for the 
motion of test particles in the resulting gravitational field as it is usually done. 

We will look for solutions of the Hamilton equations near Minkowski space-time, so that 
the matter content of the non-flat space-time is assumed to be restricted by a ultra-violet 
(UV) cutoff M: this avoids the appearance of strong gravitational fields. This UV cutoff will 
allow us to avoid the slow motion approximation and the Post-Newtonian (PN) expansions 
for the N point particles. Naturally our linearization will not be reliable if we look at 
distances di from the i — th particle of the order (or less) of the Schwarzschild gravitational 
radius Rm = 2G J /l of the particle, i.e. we must have di > Rm- Since we have a unique 
treatment for the near and far zone (we have matter everywhere) we do not need radiative 
multipoles for the gravitational field but only matter multipoles.We will do a multipolar 
expansion of the energy-momentum tensor in terms of relativistic Dixon multipoles on the 
3-space S r and we will recover the standard quadrupole emission formula. 

An important aspect of the York canonical basis is the separation of the 12 ADM equa- 
tions for the gravitational field in three sets: 

a) the four contracted Bianchi identities for the time derivatives of the unknowns in the 
super-Hamiltonian and super- momentum constraints (they imply that, if the constraints 
are solved on the Cauchy surface S To , then the solution is respected on the subsequent 
3-spaces E r>To ); 

b) the four equations for the time derivatives of the primary gauge variables (the 3- 
coordinates and the York time): since we have fixed these gauge variables these equations 
becomes equations of elliptic type for the lapse and shift functions (the secondary gauge 
variables 1 + n(r, a) and fj( )(r, <?)) of our family of gauges (instead in the harmonic gauges 
the lapse and shift functions obey wave equations); 

c) the Hamilton equations for the tidal variables R a , U a , a — 1, 2, which become hyperbolic 
equations (dP R a (j, a) = ....) for R a after the elimination of the momenta U a by inverting 
the first half of the Hamilton equations (those for d T R a ). 

To these equations we must add the super-Hamiltonian and super-momentum constraints, 
which are equations of elliptic type for their unknowns on the 3-space E r , namely for the 
conformal factor </>(r, a) of the 3-metric (namely the 3- volume element) and the off-diagonal 
terms cr^ a )(b)\a^b(T,^) of the shear of the congruence of Eulerian observers associated with 
the 3+1 splitting. 
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All the previous equations and also the matter Hamilton equations are considered on the 
instantaneous 3-spaces of the whole space-time without making a separation between the 
near and far or radiation zones with respect to the particles for the study of gravitational 
radiation. 

As a consequence, only the equations for d\ i?a(r, a) will lead to linearized equations in- 
volving the flat d'Alambertian and implying matter-dependent retarded solutions R^ 6 ^ (t, a) 
(with the homogeneous solutions eliminated by a no-incoming radiation condition with re- 
spect to the asymptotic Minkowski 4- metric). We will see that the linearized wave equation 
implies that only the TT (traceless-transverse) part of the matter stress tensor is relevant. 
Therefore we recover TT gravitational waves also in this family of non-harmonic gauges. 
Moreover, by means of a transformation of the 3-coordinates o r on the 3-space E T , we can 
find a generalized TT gauge with the relativistic inertial effects connected with the York 
time 3 K explicitly shown. 

Moreover we find the correct energy balance (the back-reaction problem) for the emission 
of GW's by using the conserved ADM energy and avoiding singular quantities like the 
gravitational self-forces on the particles [16, 17] due to the Grassmann regularization and 
objects like the Landau-Lifschitz energy-momentum pseudo-tensor of the gravitational field. 

The solutions of the linearized super-Hamiltonian and super-momentum constraints in- 
volve only instantaneous quantities on the 3-space S T and the same happens for the lin- 
earized equations for the lapse and shift functions. These are the instantaneous action- 
at-a-distance effects. Moreover the linearized Bianchi identities turn out to be automat- 
ically satisfied. Therefore, the linearized solutions for 0, cr( a )(b)|a^b, 1 + n , "-(a), depend 
upon matter-dependent instantaneous quantities and upon the instantaneous tidal variables 
= lar Ra- When we replace with the retarded solution vi 1 ^ ret ^ = 7^ r R^ et \ 
the previous variables become functions of suitable combinations of instantaneous and re- 
tarded terms all depending only on the matter. At the lowest order in 1/c 2 the results of 
the harmonic gauge used in the IAU conventions for the solar system are recovered (if 3 K 
is negligible inside the solar system). 

At this stage the extrinsic curvature tensor 3 K rs {r, a) is completely determined except 
for the numerical function describing its trace 3 K(t, a). Therefore, the final determination 
of the dynamical instantaneous 3-spaces S T (whose inner 3-curvature is determined by the 
gravitational waves and matter) associated to the linearized solution requires a choice of 
the York time: this has to be done by relying on the observational conventions, which hide 
a notion of clock synchronization. Instead in the harmonic gauges one uses the Euclidean 
instantaneous 3-spaces of the inertial observers of Minkowski space-time, avoiding to look 
at the extrinsic curvature associated to the harmonic solutions! 

Both in the previous solutions and in the resulting Hamilton equations for matter the 
influence of the relativistic inertial effects connected with the York time 3 K are explicitly 
shown. It turns out that at the HPM level all the equations depend on the following spatially 

non-local function of the York time: 3 /C d = ^ 3 K. 

In a third paper [18] we will give the PN expansion at all orders of the matter equations 
of motion (in absence of the electro-magnetic field) in the slow motion limit resulting from 
the HPM linearization. There we will evaluate the dependence on the inertial gauge variable 
York time 3 K^ = A 3 /C(i) of geometrical and physical quantities. The relevance for astro- 
physics and cosmology of these results will be discussed in the third paper, after an analysis 



6 



of the relation of the gauge problem in general relativity with the conventions used for the 
description of matter (extended bodies) in the geocentric (GCRS), barycentric (BCRS( and 
celestial (ICRS) reference frames. There we will discuss the possibility that dark matter can 
be simulated with these relativistic inertial effects in a Post-Minkowskian (PM) extension 
of ICRS. 

In Section II we review the Hamilton equations of the gravitational field and of matter 
in our family of 3-orthogonal gauges, which were given in Appendix C of paper I in the 
3-orthogonal gauges defined in that paper. 

In Section III we define our linearization scheme and we define the HPM expansion. Also 
the coordinate transformation connecting the 3-orthogonal gauges to the harmonic ones at 
the lowest HPM order is defined. 

In Section IV we solve the linearized constraints and the equations for the lapse and 
shift functions in our family of 3-orthogonal gauges. Also the ADM Poincare' generators are 
given till the second order. 

In Section V we give the linearization of the equations of motion for the particles and the 
electro-magnetic field. 

In Section VI we give the linearized second order equations of motion for the tidal vari- 
ables R a and we show that they imply the wave equation (with respect to the asymptotic 
Minkowski 4-metric) for the TT part 4 hT^ rs of the diagonal 3-metric A g(i) rs on the 3-space 
E T . Also a generalized (non 3-orthogonal) TT gauge is identified 

In Section VII we study the retarded solution for HPM GW with asymptotic background 
and we show that the dominant term is the emission quadrupole formula by using a multi- 
polar expansion of the energy-momentum tensor in terms of Dixon multipoles. Also the far 
field behavior is considered. Then we study the energy balance associated to the emission 
of GW's. Finally we look at the problem of detection of GW's. 

In the Conclusions, after a review of the results and of the problems which will appear at 
the second HPM order, we discuss the gauge problem in general relativity, the dependence 
upon the observational conventions for the 4-coordinates of the description of matter and 
the relevance of the York time 3 K^- ) for explaining at least part of dark matter as relativistic 
inertial effects (to be discussed in the third paper [18]). 

For a comparison with our formulation, in Appendix A there is a review of the standard 
approach to GW's by using Einstein's equations in harmonic gauges. 

In Appendix B there is a discussion of Dixon multipoles and of the multipolar expansion 
of the energy-momentum tensor. 

In Appendix C there is the study of the balance for momentum and angular momentum 
in GW's emission. 
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II. THE HAMILTON EQUATIONS IN THE 3-ORTHOGONAL SCHWINGER 
TIME GAUGES 

As shown in paper I, the 3-orthogonal gauges of ADM tetrad gravity are the family of 
Schwinger time gauges where we have 

a(a)(T,a)mQ, ip( a ) (t, a) « 0, 

e\r,a) «0, ^(r,3) = J^ 3 K(r,<?) k ^F^a), (2.1) 
in the York canonical basis (a, r = 1, 2, 3, a = 1, 2) 





n n (a) 


QT 




i2al 




7T„ « 7T n(a) « 


7T<*> 




nJ 



(2.2) 



The tidal variables of the gravitational field are I^a- The primary gauge variables 
(the inertial effects connected with the choice of 3-coordinates on S T and with clock syn- 
chronization) are 9 r and n^. their conjugate variables nf ^ and are the unknowns in the 
super-momentum and super-Hamiltonian constraints, respectively. While 0(r, a) is the con- 
formal factor of the 3-metric (namely the 3- volume element), tt^ may be replaced with the 
off-diagonal terms cr( a )(b) \a^b(j, <?) of the shear of the congruence of Eulerian observers associ- 
ated with the 3+1 splitting (the diagonal elements are connected to the tidal momenta n s ). 
The secondary gauge (inertial) variables (determined after a gauge fixation of the primary 
ones) are the shift {n r = 3 e^ n^) and lapse (1 + n) functions. 

The canonical variables for the particles are t)1(t), K, ir (r), i = 1, .., TV (we use the notation 
K ir instead of the one k ir of paper I). For the electromagnetic field in the radiation gauge the 
canonical variables are A± r (r, a), 7t^(t, a) and from paper I we have the following expression 
for the electro-magnetic fields F rs = d r A± s — d s A± r , B r = e ruv d u A± v , E r = —F TT = 
—d T A± r + d r A T with A T given in Eq.(3.32) of I. 

The members of this family of gauges differ for the value of the York time 3 K(t, a) ~ 
F(t, a), where F(r, a) is a numerical function (for F(r, a) = const, we get the constant mean 
curvature (CMC) gauges of ADM theory). This gauge variable, describing a relativistic 
inertial effect (which does not exist in Newton theory in Galileo space-time), is the remnant 
of the special relativistic gauge freedom in choosing the convention for clock synchronization, 
i.e. for the identification of the instantaneous 3-spaces E r . While in special relativity the 
whole extrinsic curvature tensor 3 K rs (r, a) of S T is pure gauge, in general relativity it is 
determined by the dynamics with the exception of its trace 3 . 

In this Section we will give the restriction to this family of gauges of the Hamilton 
equations and of the constraints given in paper I by using the results of its Appendix C. 



3 Let us remark that positive (negative) extrinsic curvature implies that the instantaneous 3-space, as an 
embedded 3-manifold, is a concave (convex) 3-surface of the space-time 
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We will use the following notational conventions (most of them are defined in paper I): 

a) We write the conformal factor of the 3-metric in the form 0(r, a) = 4> 1 ^{t, a) = e q( - T ' a \ 
Then we have 4>~ l d r <fi = | _1 d r 4> = d r q, _1 d 2 <fi = d 2 q + (d r q) 2 (this notation was not 
used in paper I). 

b) We use V ra for V ra (8 n ) to simplify the notation. We use the notation V ra (0) = 5 ra , 

V(i)ra = dVr QQi ' |e*=0) B(j)jw — ^—jfjjr ~ l.«=o- As sa id in Subsection IIC of paper I, we use 
angles # l corresponding to canonical coordinates of first kind on the group manifold of SO (3), 
because this implies V^ rs = 2B^ rs = e irs . 

c) The set of numerical parameters 7 Sa satisfies [6, 7] J2 U 7s« = 0, J2 U 7a«75« = 
___a 7»« lav — S uv — |. A different York canonical basis is associated to each solution of these 
equations. 

As shown in paper I, in the York canonical basis we have the following building blocks for 
the Einstein-Maxwell- Particle system in the radiation gauge of the electro-magnetic field 4 
( 3 e( a ) r and 3 e^ a - ) are cotriads and triads on S r respectively) 



4 

9tt 



4 

9rr 



= e[(l + n) 2 -$> 2 a) 



= -e ^ ^(a) 3 e( a )r = -e 1/3 Q r n (r) , 

a 

_Vs = -e 3 £„ = -e ^2 3 e(a)r 3 e( a )s = -e(j) 43 g rs = -e<t) 2/3 Q 2 r 5 rs , 



1,2 

.(i) 



3 e( a )r = 3 e( a )r = ft 1 * Qa$ra, 3 e( a) = ^ e\ a) = 1/3 Q a 1 5 ra , 



Jo )= c 3 



87TG 



^ Qa^fo 1 e ^°"(«)(&)U^&' 



aft 



87rG e a bi (e) 



c 3 1 7 8tt G 



4 The so-called gothic inverse 4-mctric h AB = \/— 4 g A g AB — ir q AB takes the form: h TT = e(j+^ — 1), 
h rr = -e Q^ 1 h(r), h rs = -e [f^ Q' 1 Qj 1 ((1 + n) 2 S rs - n {r) n {s) ) - S rs }. At spatial infinity we get 
h AB -> 0. 
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»~ r 1/3 (Q 2 r 5 rs [2 £ T6r H 5 - ^] - 2 Q P Q s - — ) = 



b 

2/3 <5r + (1 ~ <5 rs ) <7( r )( s ) Qr ~ _1 ^ 7ar ^ Q 2 . <5 r 



2 de/ 1 1 /^ 87rG A 2 I-2 |~V^ TT2 , 

• (2-3) 



a 6 



fvrf) 2 fvrf) 2 fvrf^ 

~T~ ^ 



(Q 2 Qg 1 - Q 3 Q2 1 ) 2 (Q 3 Qr 1 " Qi Qs 1 ) 2 (Qi Q2 1 - Q2 Qi 1 ) 2 

Eqs.(2.18) and (6.2) of I have been used. 9 = —e 3 K and (T( a )(b) are the expansion and 
the shear of the congruence of Eulerian observers of S r as shown in paper I. In paper I it 

is also shown that the original momenta conjugate to the cotriads 3 e^ a y before going to the 

_ ( $ ) 
York canonical basis are 3 ix\ a) « <p~ 1 ' 3 [5 ra Q" 1 (0tt^ + Y\ Taa n a ) + £• Q~ x Qr ^\ Q -^ \ 

From Eqs.(6.4) of I we have the following expressions of the mass and momentum densities 
(we use the notation M. and n ir instead of M. and k ir of I for the quantities in the electro- 
magnetic radiation gauge; see the Introduction for c(a, a )) 

M(r,a)\ ei=0 = Vi(r)) Vi ^mj c 2 + ^ £ Q~ 2 (^(t) - ^ A ±a ) * (r, a) + 

+ Yc I^ 1/3 ( ^ Ql 6ra 6sa nl nl + 1 ^ Q ~ a 2 g& ~ 2 Fab Fab ) ] (r ' ^ " 
- Yc [^ VS E « *• *- ( 2 ^"E s " m E 0- * ^|Jr^) 

rsan m j 

5 ^ — ^ — ( r ' (T )' 



(9 a™ 

3 



(2.4) 



From Eqs. (C8), (C12) and (C20) of I we get the following expressions for the Laplace- 
Beltrami operator and for the intrinsic 3-curvature 3 R[9 n , <fi, = 0~ 5 ^ — 8 A<p + 3 R<pj = 

0-6 ^ + T j with 3fl[0n = 0-2 ^ + ^ j 
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A\ ei=0 = E Qa 2 (r, a) [d 2 a - 2 ^ 75a d a R- b (r, a) d a 

a I 

= E^'fr*) \dl-2d a Y^\r,a)d a 



T(t, a) 



^=0 



= -2 



71 (r, a) 



e*=o 



3 R(r,a)\ e , =0 



1/3 E Q« 2 ( E [ E ( 2 ^ 7 C -a - <y 9 a ^ s a ^ 

a (,5 

4 75a ,9 a g9 aJ R 5 ] + 8(<9 a g) 2 ), 

1/3 E^ 2 ( 4 [^9 + 2(9«?) S 

a 

9 a 2 If> - 2 (5 «9 a q - «9 a r«) «9 a r«)] (r, a), 

-2 [0 1/3 E Q« 2 ( 4 ^) 2 - d l r i 1} - 2 (d a q - d a rf)) d a lf>)] (r, a), 

a 

E (<£ 2 [ 2 52 r i 1} - E ^) 2 - 2 r i x) ) 1 ) ^ *)■ 

a 6 

(2.5) 



In the 3-orthogonal gauges we have the following form of the constraints, of the Bianchi 
identities and of the equations of motion of gauge and physical variables. 



A. The Constraints 



The super-Hamiltonian and super-momentum constraints determining = 1 / 6 = e q and 



a (a)(b)\a 



SwG 7-1 



Q b Q?-Q a Q- b ^ giv6n in EqS - (6 - 5) ' (6 - 6) ' (C12) ' (C20) ° f are 



[Eqs.(2.4) and (2.5) are also needed; we use the notation H( a ) instead of %{a) of I] 



+ 



2%G 
8n 2 G 2 



0(r,a) ((A| ei=0 - ^ 3 R\e>=o) <P + ^r<P 1 M\ 9 i =Q + 

^ 7 E n ^ + ^ 5 E ^(.)W^(«)W-^^ 5 W) (T '^ W0 ' 



ab,a=ib 



(2.6) 
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"H(a)|^=o( r , O 5 ) 



b^a 

„ , [2 „ , r/ SttG ~ i 



2/3 (r, a) ( £ QjT 1 <7 (a)(b) + (6 d b q + d b (lf> - lf>)) <r (a) 

( £ ( 75a 9 a U- b - 9 a ^ n 5 ) + M {a) )] ) (r, a) » 0. 

(2.7) 



All the constraints depend on the York time 3 K. 



B. The Contracted Bianchi Identities 



The Hamilton equations for the unknowns in the constraints, given in Eqs.(6.7), (6.8), 
(C21), (Cll), (C19) and (C5) of I, are (= means evaluated by means of the equations of 
motion) 

d T 4>(r,a)\ ei=0 ^6(4>d T q)(r,a)\ ei=0 = | - (1 + n) 4> 3 K + 2 / 3 ]T Q' 1 (d a n {a) + 



+ n (a) (49 a g-9 a r«))](r )( 7), 



(2.8) 



d T vrf ) (r, a)\ ei=0 = e ia6 Q a Q b 1 <9 T <t (o )(6) + (6 d T q + d T (T^ - r^)) <7 (a 

ab 



*)(&) 

(2.9i 



We can extract d T a^ b ) from Eqs.(2.9) due to Eqs.(2.3), i.e. 
</> Eab CiabQaQb 1 CT( a )( b ). By using V (i)rs = 2 S (i )„ = e irs and Eqs. (6.8) and (C5), (Cll), 
(C19), (C21) of I, we get 



d T cr{a)(b)\a^b,ei=o = - [{g d T q + 



QbQa 1 + QaQ b 1 



QbQ- a l -QaQ b 

^-abi 



^ Q b Q~ l -Q a Q' 1,y " 



J d 3 a, ((l + nCr,^)) 



-° 2 <W(r,<?) lfl,=0 . 



- 
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f ja ( A , / -* \\ 5M{t,B{) 
J da, (l + nCr,^)) l^o = 

^ 2/3 E™ Qa 2 (^)ra <* sa + S ra V (i)sa ) (k^t) - f A ±r ) (^ s (t) - f A ±s ) 

+ (l + n(r, a)) (0- ^ i [ Ql V H)ra Ssa^l < + E V(i)ra F rb F ab - 

ars abr 

-\Y.Ql + ( 2 «± - E 5rm E & »* 9 ™ ^( r ))) 

arsn m i 

6Sn E QjVjdnc(a,Vj(T)) )(t,5% 



y rfV [l + n(r,5i)] 



£«S(t, gi) 



=o — 



-2 



1/3 g a 2 V (i)ra (d a n d r q - Tba d r R b ) + d r nd a q r br d a R b ) 

ra b b 

- (1 + n) [2 (2d a qd r q- d a d r q) + ^ (Tba + Tbr) d a d r R b + 

b 

+ 2 E ( Tba da Q d r R b + Tbr dr q O a R b ) - 
b 

- ^2( 2 TbrTa + 6 b5 )d a R b d r R 5 J (r,a), 



be 



d 3 ai n(r, (Ti) 



=0 — 



= -2 



1/3 g a 2 V (i)m (<9 r <9 a n - 3 (d r n d a q + d a n d r g)) (r, a) , 



13 



8ttG f 3 ^ _ | 



- !((i+»)££ 



(m)(n) 



-8ti G 



- ^ ^ tmnttikttkcdQcQd 1 °"(c)(d)) ( r , <?) 



[t~ 1/3 E ( E * *w [^r ^ ^ E (ts. - 7 5 .) n 5 



~ ^ (E Qb 1 e irb&(a)(b) + Qa 1 E ^ 1 e iafe 0"(r)(b)) 



+ 



+ 



i:i , 7 M 3 ^-E^ jR 5n 5 ) + 

a b b 

]T 7 - bb 9 6 n- 6 + (75a - Tbb) (2d bq + d b rf>) n 5 + M 6 ] + 



+ 



o- (6)( C ) - e icb d a a ( 6 )( c) + 4 (e ia6 <9 C q 

b,c^b 

e ic b d a q) <7(6) (c ) - (e iafe <9 C - e ic6 9 a ) (r^ - r[ 1} + T^) a {b){c) ^j - 
E Ql l °{a)(b)tirb{2d r q + d r T a 1] ) ) (r,a). 



C. The Shift Functions 



(2.10) 



From the r-preservation of the gauge fixings 1 (t, a) ~ 0, see Eqs.(6.10) of I, we get the 
following equations for the shift functions 



Q b 1 d b n {a) +Q a l d a n {b) - Q b 1 ^2 d b q + d b n {a) 
+ Q- 1 (2 9 a g + 9 a r[ 1) )n (6) ])(r, ( 7)^ 
« 2 1/3 (1 +n)cr (a)(6) | a ^ (t,<t), a ^6. 



(2.11) 



D. The Instantaneous 3-Space and the Lapse Functions 



The preservation in r of the gauge fixing constraint 3 K(r, a) ~ F(t, a) given in Eqs. (2.1) 
gives Eq.(6.12) of paper I. The restriction of this equation to our family of gauges (by using 
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Eqs. (C22), (CIO) and (C17) of I and the super-Hamiltonian constraint) gives the following 
equation for the determination of the lapse function (it is the Raychaudhuri equation) 



(1 + n) [$-V*M J d'a, (l +n(r,a 1 )) 5 M ^ 



a 



5 (f>(r, a) 

\ 2/3 w + ( ^) 2 r 4/3 E n ^ + ^ E 

a ab,aj^b 



- (1 + ra 

- 2 / 3 [ - 9 r 3 K + 0- 1 / 3 E *(«) 1 ^ ] ) (T, 3) = 0, 



= -2 E *V, ((1 + n) V -i -, r } ( r. a ) 



'm 2 C 2 + 0-2/3 Ea Q-2 ^. a(r ) 
(l + n(r, <?)) (> 1/2 [± E ^ W + ^ E 2 ^ F ^ " 

ars a& 

i E Ql 6 ra $sa (^1 ~ Yl ^ E ^ ^ 9 ™ C ^ 

arsn m i 

^E^^^c(^^(r))])(r,a). (2.12) 



We can also consider gauges in which in Eq.(2.1) we have 3 K (r, a) ~ F(r,a,a — 
ffi(r),...,a - ff N (r)). In this case we have d T 3 K = d T F\^ T ) + J2i f| 5 7j( r ) with ^l 7 ") 
given by the first set of Hamilton equations for the particles. 

E. The Dirac Multipliers 

Once the lapse and shift functions are known, the Dirac multipliers appearing in the 
Dirac Hamiltonian, given in Eq.(3.48) of paper I, are determined by the following equations 



d T n(r,a) = X n {r,a), 
d T n (a) (T,a) = Afi (o) (T,d% 



(2.13) 



F. The Weak ADM Energy 

The weak ADM energy, given in Eqs. (6.3) and (C8) of I, becomes (S is given in Eq.(2.5)) 
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EadmIb^o — c J d 3 cr |yW |g» =0 — — S\ 6 i =0 + 



AttG ~ i \ -v 2 

+ — ^- 1 > n 2 



a^b 



(2.14) 



where Eqs.(2.4) and (2.5) have to be used. As noted in paper I the gauge momentum 
proportional to the inertial York time 3 K (existing due to the Lorentz signature of space- 
time) gives rise to a negative gauge kinetic term, without analogue in ordinary gauge theories 
(electro-magnetism and Yang-Mills theory). 

G. The Equations of Motion for the Tidal Variables 

By using Eqs.(6.13), (C9), (C15), (C23) and (C7) of I, we get the expression of the 
momenta n a and then the second order equations of motion for the tidal variables R & - For 
the tidal momenta we get 



n a (r,£) 



8nG 



8ttG 



7a 



a 

1 + n r, a) L *-^> 

(2d a q + d a r^) - d a R^ U(a) - 7aa da 71(a)) (t, 0% 



(2.15) 



As a consequence, for the tidal variables R a we get (Eq.(2.15) is used to eliminate the 
dependence upon II a of the equation of paper I) 



d 2 T R- a (T, a)= « [4>- 1/3 n { a) £ ( 7 aa 75a " ^) 5a 9 T i2g + 

a 5 

+ 0~ 1/3 ]T Qa" 1 [(Too (2 <9 a g + 9 a r«) - d a R-^j n {a) - 

a 

~ 7aa<9 a n (a) ] 9 r rW _ 

- 0' 1 i? a + \ 0- 1/3 £ Qa" 1 ( [ina (~d a q + d a T«) 

a 

- <9 a i? a ] n (a ) - 7 aa <9 a W( a )) ] <9 T 0- 
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+ 



a 

d T R- a + 0- 1 / 3 J2 Qa 1 ( [Taa (2 d a q + d a if } ) 



- d a R- a \ n (o ) - 7aa <9 a n (a ) 



<9 T n 
1 + n 

- F 1/3 E 1 ( [t*« (2 <9 a g + <9 a lf>) - d a R- a 

a 

- laad a d T n {a) ) (t,<t)+ 



d T ri(a) 



+ ^(^(l + n))^) fd'a, [(l + n)(r, 



^ <J5(r,5i). , . JT(r,?i; 

^l) 7TT/ =h\8*=0 + n \ T i 



5i?a(r, a) 5Ra(r,a 
/0 2/3 r/„ _ _ ... . „™ 0„n 



le 4 =o 



+ 



+ (1 + n)) (r, a) E Qa" 1 [(9. n ( «) + *(.) (4 9 a g - 9 a T« - ^J-)) (d r Jfe + 

a 

+ ^ 1/3 Eg; 1 [(7. c (2«9 c g + ,9 c r( 1 )) -9 c i? s )n (c) - 7sc 9 c n (c) ]) + 

c 

+ n (o ) <9 a (<9 r R n + 0~ 1/3 53 Q7 1 (jac (2 <9 C g + d c T^) - d c R^j n (c) - 7sc d c n (c) 

c 

2/3 E ilaa - lab) [d b W(a) ~(2d b q + d b T^) n (a ) (T(a)(fe)) (t, a) - 



+ 



|6» l =0 



= 2 (0 1 / 3 J] g- 2 [d a n (2 7sa 9 a 9 - J] (2 7sa 7ba - fo) d a - 

a I 

- (1 + n) (2 7 , a (-,9 2 g + 2( ( 9 a g) 2 ) + 

+ E( 2 7«« 75a - (<9 2 ^5 + 2 <9 a q d a R~ b ) + 

b 

+ 53 ( 2 Tba ^ ~ iaa 7ea ) ~ %a ^) 5a R ~ b da Re ) ) ( r > ^ ' 



6c 



/ 



rf <7in(r, cti) — — lfli=n = 2 



<5-R a (r, it) 



1/3 53 ^ n-Gd a qd a n)] (r, a) , 
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y rf 3 <7i ( 1+ „ (Ti(7i) ^_^_i = 

= -^f'(M,(r))»((l + ») 

i 

^ m 2 C 2 + 0-2/3 Q-2 ^. a(r) - & 

+ (l + n(r, <?)) [I £ 7sa ^ ^ ^ ^ _ J_ ^ (%a + 7sfc) Q -2 q-2 Kb F( 

ars ab 

-^7« <*ra 5sa (2 Til - E ^ E ^ ^ 9 ™ ^W)) 

arsn rrt i 



(2.16) 

The three integrals at the end of Eq.(2.16) were given in Eqs.(C28), (C29) and (C27) of 
paper I. The expression of the last integral in Eqs.(2.16) has been obtained by using the 
super- momentum constraints (2.7). 

To get the final form of the second order equations for R a we have to use: a) Eq.(2.8) for 
d T 4>; b) the Hamilton equations (2.13) for the Dirac multipliers in the 3-orthogonal gauges 
with d T n and d T n^ r ) determined by the solution of Eqs. (2.12) and (2.11); c) Eq.(2.15) for 



n a . 



H. The Equations of Motion for the Particles 

By using Eqs. (6. 14) and (6.15) of paper I, the first half of the Hamilton equations for the 
particles implies 

•TV \ / V 



J m 2 C 2 + 0-2/3 Ec Q-2 ^ jc(r) _ S A±c )j 2 

Ktr (r) = ^A ±r (T,rf t (r))+m l c ^f) 2/3 Q 2 r ^[(r) + 0~ 1/3 Q; 1 n (r) ) ((} + n)" - 
/ ~ \ 2 \ —1/2-1 

E^(^ c(r)+r 1/3 "w) ) ]( r >^ r ))- ( 2 - 17 ) 



c 



18 



so that the second half of the Hamilton equations becomes 



d 



(rriiC ( 



n(r) 



^(l+n) 2 - 0V3 Ec Q2 ^c (r) " ^-1/3 Q -l ^ 

) +^^r)(r,rfi(r)), 



)(r,^(r)) = 



(9 ?? f r ' ' c vn v ' ' d dr 



W = 



/ [(1 + n) W (n) + f 1/3 E Q: 1 n w W„] (r, a), 

a V i 



5™ E Qi n 



dc{a,ffj{T)) 



W r (r,a) = -- Frs(r,a)5 sn E Q 



iVi 



dc(a,ffi(T)) 



1 + n) - ~^ J2 Ql (#00 + r 1/3 Q' c l n {c) ) 

c 

_ (1 + n) |!L + ^ EQ J^l_ 

- (2d r q + J2 Taa d r R a ) n (a) ) (t)?(t) + 0- 1 / 3 g; 1 n (o) ) + 



-1/2 



+ 2/3 E G« ( 2 9 ^ + E ^ 9 - ^) (A a ( r ) + ^ V3 G« 1 "(«: 



(2.18) 

Here W is the non-inertial Coulomb potential, F ir are inertial relativistic forces and the 
other terms correspond to the non-inertial Lorentz force [4]. 

I. The Equations of Motion for the Transverse Electro-Magnetic Field 



Finally, from Eqs.(6.16) of I the Hamilton equations for the transverse electro-magnetic 
fields in the radiation gauge become (P±(a) = 5 rs - J2 U v ^™ $ sv ^r) 
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-1/3 



d T A ±r (r,a) = ^^^(a 

nua 

m a i 

r 2/3 (i + ») g„- 2 ««,m 



+ «- («1 - E «" E «- * 8c( a^ (r)) ) + 



[ 



m? c 2 + 0- 2 /3 £ 6 g- 2 (^ 6 (r) - f A ±b )' 



2 (i + n ) Qf Qf 5 ma (d b F ab - [2 d b q + 2 d b (r« + if*) 

+ 2 0- 1 / 3 ^ g; 2 g 6 ~ 2 5 ma d b n F ab - 



F ab ) + 



ab 



1/3 E "w Qa 1 < - [2 9. « + <9 a r« 



^ m 1 

7T, + 



1 

7T, + 



+ 6 ma 

n 



E (K< 



dc(a,ffi(T))) _ <9 2 c(a,^(r) ))- 
9cr n da n da n 

dc(a,fh(r)) 



+ 



+ <t> 1 l^Qa (d a n {a) — 5 ma ^ d n n (a) — 



(2.19) 



J. The Weak ADM Poincare' Generators 



While the weak ADM energy P^dm ~ \ Eadm is given in Eq.(2.14), Eqs.(2.22) and (3.47) 
of paper I give the following expressions in the 3-orthogonal gauges for the other weak ADM 
Poincare' generators (the last term in the boosts was added in Ref.[19]) 
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ADM 



= J d 3 a [ 3 g rs M s - 2 3 r su (r, a) 3 IT«] (r, a) = 

= 2 I d 3 a { 0- 2 / 3 Qr 2 ( 2 75, d r q + £(7Sr7*r - \ hi) ft-ife) U~ b 

J b a - 

- T ^^ /3 Q; 2 ^d rq + d r r^) 3 K + 

J2 Q^Qt (2d dq + d d if>) (7 (r)(<0 + \r 2/3 q; 2 m t } 



87TG 



jrs 

J ADM 



= 2 fd 3 a{ a r [ 0- 2 / 3 ( 2 7 5s d s q + ^( 75s7 , s - (1/2)%) d s R- a 

~ T ^^ /3 Qf(^ s q + d s T^) 3 K + 



+ 



8nG 



- a- 



+ 



0- 2/3 £ ( 2 7" br 9 r g + Y^iTbrTar - (1/2)%) 9 rJ R s ) n- b - 

b a 

-^4> 1/3 Q; 2 (4d r q + d r r^) 3 K + 

£ Qr'QJ 1 ( 2 ^ ? + a* if>) (7 (r)(<0 + ^- 2/3 g r - 2 At ] }, 



8ttG 



jtt jrr 

J ADM — J ADM 

8tt G 



d 3 a a r 



16nG 



r~ 3 ns /3-pu 3-pv 3-pu 3-pv \ 

V 7 9 \ nv su ns vu) 



+ 



c 3 v/7 

c 3 
16nG 



'G nsuv 3 Yl ns3 U uv -M] + 
51 ( 3 g vs - 5 VS ) d n [^7 ( 3 9 ns 3 9 uv ~ 3 9™ 3 9™)] ) (r, a) 
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/''Mub*-*-^- 1 ^- 

b 



167rG r 

s 

- «9 r (r« + r« + g)] }(r,*)«0. 



(2.20) 



As discussed in Section HE of paper I, JYdm ~ are the gauge-fixings for the rest-frame 
constraints P^dm ~ eliminating the internal 3-center of mass in the 3-spaces S r , which 
are non-inertial rest frames of the 3-universe. 



K. Dimensions 



In checking the validity of the previous formulas it is useful to remember the dimensions 
of the relevant quantities: 

[ T = d] = M = [a] = m = [I], m = [m, t c] = [P»] = [E/c] = [mlt' 1 }, [S] = [h] = 
[J**] = [mlH-\ [T AB ] = [M] = [M r ] = [H] = [H {a )\ = [ml-H-\ [*g] = [*g] = [n] = 
[«(«)] = [\a)r] = m = m = [$] = [0], [ 3 ^ a) ] = [ 3 n-] = [U- a ] = [^] = [ml-H~\ 
[V (a) ] = [ 3 ^ rs ] = [ 3 K] = [a {am ] = [I' 1 ], [ 3 R] = [ 3 Q rs{a) ] = [l-% [Q,] = [m^r^t-% 
[A ±r ] = [Ql- 1 ] = [m 1 ' 2 VH- 1 ], [*A ±r ] = [mlt- 1 ], K] = [Eft = [B r ] = [l' 1 A ±r ] = 
j m i/2 r i/2 r i] > [G = Q.710- 8 cm 3 s- 2 g- 1 ] = [m^lH' 2 ], [G/c 3 = 2.5 1(T 39 sec/g] = [m' 1 1], 
[G/c 2 = 7.421 1(T 29 cm/g] = [rrT 1 1]. 
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III. THE WEAK FIELD APPROXIMATION AND THE LINEARIZATION 



The standard decomposition used for the weak field approximation in the harmonic gauges 

is 

9iu> Vfiv h^iv ) I h^i/ 1 , 1 a h^ v | , 1 a Op h^ v | < < 1 , 

(3.1) 

where A n lJiV is the flat metric in an inertial frame of the background Minkowski space-time. 
This is equivalent to take a 3+1 splitting of our space-time with an inertial foliation, having 
Euclidean instantaneous 3-spaces, against the equivalence principle and against the fact 
(explicitly shown in paper I) that each solution of Einstein's equations has an associated 
dynamically selected preferred 3+1 splitting. 

In this Section we shall define a linearization of Hamilton-Dirac equations in the (non- 
harmonic) 3-orthogonal Schwinger time gauges (2.1) using as background the asymptotic 
Minkowski 4-metric existing in our asymptotically Minkowskian space-times. Actually we 
look for the following decomposition to be done by using radar 4-coordinates adapted to an 
admissible (see Ref.[4]) 3+1 splitting of space-time 

A 9AB(r,a r ) = A g { i )AB (T,a r ) + 0(C 2 ) -»■ A ilAB(as ym ) at spatial infinity , 

9{l)AB(T , O" 7 *) = 4 VAB(asym) + A h(i) AB {j , a"'), 

A h(i) AB (r, a r ) = 0(() — > at spatial infinity, (3.2) 

where << 1 is a small a-dimensional parameter, the small perturbation A h^ A B has no 
intrinsic meaning in the bulk and 3 g(i) rs ( T , °" r ) = — e 4 9(i)rs( T , < jV ) is the positive-definite 3- 
metric on the instantaneous 3-space E T . In our case the instantaneous 3-spaces will deviated 
from flat Euclidean 3-spaces by curvature effects of order 0((), in accord with the equivalence 
principle. 

We must make an assumption on the variables R a , d T Ra, Ra, (ft, n, n( a ), <7( a )(b) \a^b, n s 
(or <7( )( a )), 7r^ = 12 ^ G 3 K, rjl, n ir , A± r , w r ± such that the coupled equations (2.8) and (2.9) 
(contracted Bianchi identities), (2.11) (shift determination), (2.12) (lapse determination with 
free 3 K), (2.15) (expression of Us), (2.16) (equations for d^R a ), (2.17) (expression of «j r ), 
(2.18) (equations for r}\), (2.19) (Hamilton equations for A± r and 7Tj_), all have the two 
members consistent (of the same order). 

Let us remember [14] that to avoid coordinate singularities we must always have N(r, a) = 
1 + n(r, a) > (3-spaces at different times do not intersect each other), e A g TT {r, a) > (no 
rotating disk pathology) and 3 g rs (r,a) with three distinct positive eigenvalues. 

A. A Consistent Hamiltonian Linearization 

Let us see which assumptions are needed to get Eq.(3.2) in the 3-orthogonal Schwinger 
time gauges (2.1). 



23 



The first assumption is that on each instantaneous 3-space S T we have the following 
limitation of the a-dimensional configurational tidal variables R a in the York canonical basis 



\R a (r,a) = R (1)a (T,a)\ = 0(C) « 1, 

\d u R a (r,a)\ ~ io(C), \d u d v R a (r,a)\ ~ -Lo(C), 
|9 T i2a| = -^0(C), |^i2a| = ^O(0, |9 r «9 uJ R s 
Q a (r, <t) = t-^(^) = 1 + r«(r, a) + 0(C 2 ), 

ri 1} = E ^ = °(0, E r « } = °> ^ = E ^ r ^ ( 3 - 3 ) 

a a a 

where L is a fop enough characteristic length interpretable as the reduced wavelength X/2n 
of the resulting GW's. Therefore the tidal variables R a are slowly varying over the length L 
and times L/c. This also implies that the Riemann tensor a Rabcd, the Ricci tensor a Rab 
and the scalar 4-curvature A R behave as 0(C)- Also the intrinsic 3-curvature scalar of the 
instantaneous 3-spaces £ r , given in Eqs.(2.5), is of order jp. 0(C)- To simplify the notation 
we use R a for Rm a i n the rest of the paper. 

As a consequence of the behavior of the Riemann tensor, the mean radius of curvature 
4 1Z of space-time is of order A 1Z~ 2 ps -j? 0(C). Therefore we get that the requirements of the 
weak field approximation are satisfied: 

i) A = 0(C), if A ~ R a is the amplitude of the GW; 

ii) (|) 2 = 0(C), namely L » A « 4 ft. 

As a first attempt let us put = 1 / 6 = 1 + 4>( ) + 0(i) + 0(C 2 ), n = + + 0(C 2 ) 
and n (a) = n (o)(a) +n ( i)( a ) +0(C 2 ), with (o) , n (o) , n (o) ( a ) = 0(1), (1) , n ( i), n (1)(a) ~ 0(C), and 
with similar expansions for the other variables like C( a )(&), 3 K, .... However, this implies 5 



= 7,0(0, 



,o 1 + n^ ( i) + 0(C 2 ), 



5 For |0 (1) | < 1 we have 0" = 0™/ 6 = (1 + (o) )" + n(l + 0(c))"" 1 0(1) 
0"" - r" /6 = (l + ^(o))- ,l -«(l + 0( o ))-"- 1 ^(i) ^ 0(o)= o l-n0 (1) +O(C 2 ), d r q = 0" 1 C> r = 

\ 0- 1 d r = (i + ^(o)) -1 9 r (o) + (i + 0( o ;>r 2 (i + 0( O )) 9 r <£(i) - ^(i) a r (o) i-^ 0(o)=o a,. 0(i) + o(c 2 )- 
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- e 4 g rs = 3 g rs = 4 (l + 2 ^ Tar R^j 5 rs + 0(( 2 ) = 

a 

= ((1 + (o) ) 4 + 4 (1 + (o) ) 3 (1) + 2(1 + (o) ) 4 ]T 7 , r ifc) 5 rs + 0(C 2 ) = 

a 

= (l + (o) ) 4 5 rs + O(C), 
eV T = (l+n) 2 -^4 ) = (l + n (o) ) 2 -^4 )(a) +0(C), 

a a 

-e V r = <p 2 Q r n {r) = (jl + (o) ) 2 + 2(1 + (o) ) (1) + 

+ (1 + 0( o )) 2 ^ 7 ar Ra) W(r) + O(C') = (1 + 0(o) ) 2 n (o)(r) + 0(C), (3.4) 

a 

and the equations for 0( o ), n( D ), ^( )(a) turn out to be not linear. 

Therefore we must assume 0( o ) = ri( ) = ^(o)(a) — 0. In this way Eq.(3.2) can be imple- 
mented in the following way 6 

1 + (1) +O(C 2 ), 0=1 + 60 (1) + O(C 2 ), 
1 + n = 1 + n ( i) + 0(C 2 ), n (o) = n ( i )(a) + 0(C 2 ), 

JJ- 4 9(l)AB = A VAB(asym) + 4 ^(1)AB, 

2en (1) = O(0, 
-en ( i)( r ) = 0(C), 

-2 e (r( 1 ) + 2 (1) )5 rs = o(C), if> = £ 7ar ife, ^r« = o., 

a r 

h w = e 4 r) AB4 h {1)AB = 2 (n (1) - 6 (1) ) = 0(C), (3.5) 

while the triads and cotriads become 3 £ r mt a \ — 8 r a {l — — 2 0(d) + 0(C 2 ) and 3 e(i)( a ) r = 
5 ra (1 + ri 1] + 2 0(i)) + 0(C 2 ), respectively. Therefore we have 4 g TT = e [1 + 2n {1) ] + 0(C 2 ), 
V = -en ( i)(r) + 0(C 2 ), 4 g rs = ~^9rs = ~eS rs [1 + 2 (if > + 20 (1) )] + 0(C 2 ), = 6 = 

v ^it^-;=i + 60( 1) + O(C 2 ). 

With these assumptions Eq.(2.3) implies 



6 The " trace reversed" perturbation is 4 /i(i)AB = 4 ^(i)ab — \ A f}AB e fyi), fyi) = — fyi), 4 /i(i) IT = e (6 0(i) — 

"(1)), 4 ^(l)rr = -£"(l)(r)) 4 fyl)rs = ~C [2 (if 5 +5 (1) ) -77, (1) ]<5 rs . 



= 

N = 



h(l) TT = 

h(l) Tr = 
[ h(i) rs = 
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^^n s (T,<7) = ^^U(i)a(r,a) = jO(C)=[d T R a -^2 laad a n {1)(a) (r,a) + jO(( 2 ), 

a 

8lT G 1 r^l 

V(a)(a) = <T(l)(a)(a) = 2^ ^ + £ °vC )■ (3-6) 

a 

Let us remark that everywhere appears in the combination ^§ 11(1)6 — ^ 0((), which 

behaves like <9 T i-e. it varies slowly over L. 

Finally the super-momentum constraints (2.7), Eqs.(2.3) and dimensional arguments re- 
quire 

<T(a){b)\afr = <T(l)(o)(6)|o#6 = yO(0, 



87TG 



4 e) = \ o(c 2 ) = £ (r« - if) e Mb <7 (1)(a)(6) + 1 0(C 3 ), 



3r ^ 12ttG 3 ^ 127rG 1 



4 



= 3 ^ ( i )rs = ^o(C) 



"1 111 

= (1 - 5 rs ) <7(i )(r ) W + <5 rs [- 3 K (1) - d T rW + ^ (5 ra - -) <9 a n (1)(a) J + - 0(( 2 ). 

a 

(3.7) 

These equations imply that once we have found a solution 0"(i)( a )(&) lai/j, (?) of the lineariza- 
tion of the super-momentum constraints (2.7), then we can put 7r^(r, a) ~ in the York 
canonical basis, which becomes adapted to 13 of the 14 constraints after the linearization. 

Let us remark that the triad momenta and the standard ADM momenta, given after 
Eq.(2.22) of paper I, have the following weak field limit: 7iT x = 5 ra + J2h JBa^-(i)b) ~ 

^(i-U-(i) W w+io(C 2 ), 3 n- = U%) 3 na)+%) 3 ^ = U%)+^)) + tO(C 2 ) = 



Let us now consider our matter, i.e. positive-energy scalar particles and the transverse 
electro-magnetic field in the radiation gauge. 

For the particles we have r\\ = 0(1) and f\\ = 0(1) (since r = ct, in the non-relativistic 
limit we have 7^ = Vi/c = 0(1) — 7- c ^oo 0). 

However, without further restrictions on the masses, the momenta and the electro- 
magnetic field Eqs.(2.4) would imply M = M {o) + M m + f§0(( 2 ), with M {o ) = 0(1), 
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M(i) = 0((), and M r = M.(o)r — jjO{1) (M and M r are densities; m is a typ- 
ical particle mass). But then Eqs.(2.6) and (2.7) for the super-Hamiltonian and super- 
momentum constraints would not be consistent. For instance Eq.(2.6), whose unknown is 
0=1 + + 0(C 2 ), would be M (o) (r,a) + F (1) [<j> w , M w , ...](r, a) + ff 0(( 2 ) » with 

f(i) = ^o(C). 

To get a consistent approximation we must introduce a ultraviolet cutoff M on the masses 
and momenta of the particles and on the electro-magnetic field so that 

M(r,a) = M^ V \r,a)+n {2) (T,a), 

mi = MO((), J d 3 aM%P(r,a) = McO((), J d 3 *n (2) (r, a) = McO(( 2 ), 

M r (r,a) = M {1)r (r,a), J d 3 a M$?(t,B) = McO(C). (3.8) 

Here M is a finite mass defining the ultraviolet cutoff: M c 2 gives an estimate of the weak 
ADM energy of the 3-universe contained in the instantaneous 3-spaces S T , because it can 
be assumed to be of the order of the mass Casimir of the asymptotic ADM Poincare' group. 
The associated length scale is the gravitational radius Rm = 2M ^ ps 10~ 29 M 7 . 

Therefore the description of particles in our approximation will be reliable only if their 
masses and momenta are less of McO(() and at distances r from the particles satisfying 
r > Rm (that is at each instant we must enclose each particle in a sphere of radius Rm and 
our approximation is not valid inside these spheres). This will be clear in the next Section, 
where we will obtain an equation like A^jfr, a) = j^O{Q ~ ^ M.(i)(t, a) + ... implying 
(1) (r,a) = 0(C) ~ Spafe^ + - ~ W^m°(0 namely \a-ff t (r)\ » R M . Therefore 
our results in the weak field approximation can be trusted till a distance d » Rm from the 
particles. 

From Eq.(2.4) we get for the mass density and mass current density 



7 The Earth mass M Ear th = 5.98 1 28 g gives rise to a gravitational radius REarth = 2GM j^ th _ q ; 
By comparison the Compton wavelength of an electron is — — = 3.861592 10 _11 cm, the classical electron 

radius is = 2.81794 10" 13 cm and the Planck length is L P = x f^ = 1.616 1(T 33 cto. 
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M (uv \r,a) = M^ ) V \r,a)+M^ ) V \r,a)+'}Z {3) (T,a), J d 3 aK(r,a) = McO(( 3 ) 

^([E(w)-(-i-E«.,.^™) 

a i 

^ dc{a,fij{T)) \ 1 ^ 2 l \ 



-2 (1) (^(r) - g 1 ± ) 2 - E a ri 1} (Ur) - g A ±a ) : 

E % * + E [*» + r« + if] ) (x, ,-), 



da a 

/ A /< 



<9(X a 

r/V.M^V,*) = McO(C), / dVM^V.ni = M, 0(C) 



(3.9) 



■ 1 ^ 

- i E Fr,(r,3) ( <{ r,B) - E *» E Q'* ?f ^). 

s n i 

J d 3 aM^ V r \r,a) = McO((). (3.10) 

Therefore for the particles and the transverse electro-magnetic field the validity of the 
weak field approximation requires 
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I 4 ±r (T,a , ),7r:(T,?)=0(l), with ^ -4 ±r (r, jft(r)) = McO(Q, 

c 

[ d 3 a[-7r r ± (r,a)} 2 , -F 2 s (r,a), [- F rs 7r s ± ](r, a) = McO((). (3.11) 
J c c c 



Moreover the boundary conditions at spatial infinity and the local intensities for the 
transverse electro-magnetic field must be such that the integral conditions in Eqs.(3.8) 
hold. The last line of Eqs.(3.11) agrees with Eq.(3.8) only if the radiation part of the 
transverse electro-magnetic field is concentrated in small volumes V = V 0(() with a suf- 
ficiently rapid decay outside them. The restriction — A± r (r, ^(t)) = McO((), dictated 
by Eq.(2.17), implies a bound on the value of the electric charges of the Lienard-Wiechert 
transverse potential evaluated in Ref. [15] in special relativity: since this potential has the 
form ^ A± rLW (r, a) = Ylj^i ciJ\s-^-{T)\ ^jr wrfc h F jr = 0(1), it turns out that for distances 
\fji{r) — fjj{r)\ > Rm the restriction implies that the product e^ej of the electric charges 
semiclassically simulated by the Grassmann variables QiQj must satisfy ej < RmMc 2 . 

Our results will be equivalent to a re-summation of the PN expansions valid for small 
rest masses still having relativistic velocities (^3 = 0(1), ^ = 0(1)). 

Let us remark that in this way the energy-momentum tensor T AB (its expression derives 
from Eq.(3.11) of paper I after having expressed the metric components in the York canonical 
basis, after having expressed the electro-magnetic field in the radiation gauge as said after 
Eq.(3.35) of paper I and after having done the weak field approximation) has the following 
behavior 



rpT 



M (UV) V rr 



(/t ; ir 



— A 



{^is A 



+ 



~ 7T 



m 2 C 2 + Ea(^a-^A ±a y 

dc(a,ffi(r)) 



+ 



da" 



+ 2 {2^(*±-2^ViQi — — ) 

a i 
(1) "+" /C (2)) 



Q: 



+ 



1 



ab 



29 



E T (l) = E *V - *O0) ^ y n + 

^[E«) 2 ^E^- 

a a& 

- E ^ ( 2 ^ ~ E ^ ^ 5a ^( r ))) da c (^' ^( r )) 



* — . v * — ■ v o 7ar ( K ir * ^4_Lr ) 

E^ T (i) = E 5 ^-^ r ))^ 



^E^ [w) 2 -E^- 

r a 

^ ^ Qi (2 < - E % Qj d r c (^> r /i( r ))) c(a, ffi(r)) , 

i j^i 

a i 

- E^ + l^^l- ( 3 - 12 ) 



ab 



Since, as said in Subsection HE of paper I, we have T ab {t, <r)=0 from the Bianchi 
identities and since 4 <?ab = ir nAB(asym) + 0(C), we must have 8a T^ b (t, a)=0 + 8aTZ^ ■ At 
the lowest order this implies 



8 T Mf,p + 8 r Mf^ ] = + 8 A Ufa, 



c^S + a 7ft = + (3.13) 



as in inertial frames in Minkowski space-time. The equation 8a T^(t, <j) — + 8aTZ^ 
implies 8a (t^ b (t, a) a c — Tff (t, a) a B ^j = + dA7Zf 2 ) C (angular momentum conservation). 

Finally let us consider Einstein's equations in radar 4-coordinates, i.e. a Rab — 
\ A g A B A R=^r-T A B- Since we have A R AB = t?0{C) from Eqs.(3.3) and 
^ / d 3 aT {1) A B (T,a) m R m O(() from Eqs. (3.8) and (3.12), we get from Einstein's equa- 
tions the following local estimate of the order of T^ A b(t, o) 



30 



T(i)ab(t, a) 



Mc 



0(C). 



(3.14) 



Rm L 2 



Therefore the support V = V o 0(() of the radiation part of the electro-magnetic field, 
defined after Eq.(3.11), must have V D ~ R M L 2 with L > R M . 

In conclusion, since the weak field linearized solution can be trusted only at distances 
d » Rm from the particles, the GW's described by our linearization must have a wavelength 
satisfying A ~ L > d » R M (with the weak field approximation we have A << 4 1Z without 
the slow motion assumption). 

If all the particles are contained in a compact set of radius l c (the source), the frequency 
v — j- of the emitted GW's will be of the order of the typical frequency u s of the motion 
inside the source, where the typical velocities are of the order v ~ u s l c . As a consequence 
we get v = j u s ps v/l c or A w ^ l c » Rm, so that we get - c w j << and / c >> i?M 



If the velocities of the particles become non-relativistic, i.e. in the slow motion regime 
with v « c (for binary systems with total mass m and held together by weak gravitational 

forces we have also - m J « 1), we have A >> l c and we can have l c R M - 



B. The Linearization Interpreted as the First Term of a Hamiltonian Post- 
Minkowskian Expansion in the Non-Harmonic 3-Orthogonal Gauges 

In our class of asymptotically flat space-times near Minkowski space-time the above 
linearization can be interpreted as the first term in a HPM expansion with a UV cutoff on 
the matter. This is due to the fact that most of the canonical variables in the York basis 
parametrize deviations from Minkowski space-time with the Cartesian 4-coordinates of an 
inertial frame, which vanish if G — > (for G = we get Ra = n = ni a \ — IT S = 0, = 1 in 
the Minkowski rest-frame instant form, where we have also 3 K = cr( a )(b) = 0)- 

As a consequence, by using Eqs.(2.3) we can write the following HPM expansions (here 
we do not use R a to mean R^a as is done in the weak field approximation) 




0(1). 
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Ra — ^ G n R[ n ]a, R(l)a — G R[l]a, 

n=l 

oo 

n = Yl G " n M' n (i) = 

71=1 

OO 

= E^ n "W(«)> ^(i)(a) = Cn[i]( a ), 

n=l 

oo 

= e * = 1 / 6 = 1 + G^lnh = 



n=l 

oo 



n=l 

3 00 
n=2 

Q /-I (0) °° 

^(a)(6)|a#6 = -J0 Z^TTTR n n-1 = 2^ G ^[n](o)(b)|a#b, <7(1) (a) (6) |a#6 = G <7[i] ( a )( 6 ) | 

(3 

Moreover Eq.(3.6) implies Gila = Et?=i G n H[ n ]a, so that from Eq.(2.3) we get C( a )(a) — 
|r0~ 6 Ea laaGUa = J2n=i G ™ ff N(o)(a) • Finally, since we have QaQ^ 1 - QbQ" 1 ->g^o 

(9) /flN 

2 Ea(7aa - Tab) G i2[i] S , from G <7[i] (o)(6 ) | a +b = # Ei Efl J,-^) fl [1]a We als ° § et ^ = 
V°° G™7T (e) 

The study of HPM at the second order will be done in a future paper. 



C. The HPM Linearization of the Gauge Fixings for Harmonic Gauges 

In Eqs. (5.3) and (5.4) of paper I we expressed the Hamiltonian version of the gauge- 
fixing constraints x r ( r ) <?) ~ and x r (r, a) ~ selecting the family of 4-harmonic gauges in 
the York canonical basis. 

If we eliminate the gauge fixing 9 % (r, a) ps identifying the family of 3-orthogonal gauges 
and we consider the angles (0(3) canonical coordinates of first kind) 9 z (t, a) as small quan- 
tities 0%\(t, a) = 0(() (so that we have V sa {p i ) = 5 sa + e sa j 0^ + 0(( 2 )), we can extend 
our HPM linearization to arbitrary gauges. With some calculatins it can be checked that 
the first half of Eqs. (3.7), regarding the super-momentum constraints, are still valid with 

^ = ^ 1} + o(c 2 )^o. 

Then the linearization of the harmonic gauge fixing x T (r, c?) ~ (see the second half of 
Eqs. (5.3) of paper I) becomes 

d T w(i) (t, a) « - ( Y & n(i)( r ) + 3 K (1) J (r, a) . (3.16) 
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Instead the linearization of the harmonic gauge fixings x r ( r ) (?) ~ of Eqs. (5.4) of paper 
I becomes 

d T n ( i)(r)(r, a) « - (d r n (1) + 2 d r (0 (1) - if))) (r, a). (3.17) 

Both the equations do not depend on #(i)(r, a). As said in paper I all the gauge fixings 
for the gauge variables 9^ (r, a) and 3 if(i) (r, a) compatible with these equations identify 
4-harmonic gauges. 

The previous two equations imply the following wave equations for the lapse and shift 
functions 



□ n (1) (r,a)^(2^9 r 2 (0 (1) -r( 1 )))(r, ( ?), 

r 

□ n ( i)( r )(r,(T) -d r ^2 n W(s)( T ^) ~ 

» (^drdriT^ - (1) ) + d r 3 K {1) yr,a). (3.18) 



The solution of these hyperbolic equations requires initial data at r — > —oo. Instead in 
the family of 3-orthogonal gauges we have elliptic equations on a fixed 3-space S T requiring 
data only on it. 

If we denote n^f^ and the retarded solutions of Eqs. (3. 18), the 4-metric 

4 9(i)ab(^i ^) * n harmonic radar 4-coordinates (f, a) will have the form 4 <7(i).f^ = e ^1 + 

2n \i) H) )i 4 9{?"r = ~ efl {?)(r)> ^Wrs = ~ e + A {f)fs) with A ^)rs depending on which 
harmonic gauge one chooses. The connection to the 4-metric 4 5 , (i) j 4b( t , ?) i n the family 
of 3-orthogonal gauges is by means of a 4-coordinate transformation r = f + a(i)(f,(j), 

a r = a r + b r (1) (f,a) implying *g$%b(f, a) = §§^ 4 #(i)cd(t, a). As a consequence one 
gets 



1HH) 



da {1) (T,a) . ... 
qZ = \n\ x) ~ n (1) 

9%(f,^) (_ {HH ) _ V- ^ 
= ^(l)(r; I -«(l)(r)J(T,^), 

f (f. a) db s (1) (f,a) 



df 

'(!)!_ 

9a s 



(2(r( 1 ) + 2 (1) )5 rs -A™)(f^). 



(3.19) 
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IV. THE SOLUTION OF THE LINEARIZED EQUATIONS FOR 0, 1 + n, vrf } , n (a) 
AND THE LINEARIZED ADM GENERATORS 

In this Section we find the linearization of the super-Hamiltonian and super-momentum 
constraints and then of the equations determining the lapse and shift functions of our family 
of 3-orthogonal Schwinger time gauges. The solutions of these linearized equations allow 

~ (0) 

to express 0(i), 1 + n^, <7(i)( a )(b) |a^& (i- e - ), ^(i)(o) i n terms of matter and of the tidal 
variables = ^2- ^arRa- in the next Section we will see that the linearized equations 
for the tidal variables depend only on matter, so that at the end the previous solutions 
will depend only on the matter. While the equations of elliptic type solved in this Section 
will determine the instantaneous inertial dependence on the matter of the 4-metric (like the 
Coulomb potential in the radiation gauge in the case of the electro- magnetic field), the wave 
equations for the tidal variables in the next Section will determine the retarded dependence 
on matter of the 4-metric. We will also see that at this order the contracted Bianchi identities 
are identically satisfied. Then we will evaluate the asymptotic ADM Poincare' generators 
till the second order. 

A. The Super-Hamiltonian Constraint and the Lapse Function 

Let us first consider the determination of = 1 / 6 = 1 + 0(i) + 0(( 2 ) by using the 
super-Hamiltonian constraint (2.6) and of the lapse function by means of Eq. (2.12). 

1. The Equation for <p^ 

Since the linearized Laplace-Beltrami operator and the scalar 3-curvature in the family of 
3-orthogonal gauges, see Eq.(2.5), have the expressions (A = J2 a d\ is the flat asymptotic 
Laplacian) 

A| fli=0 = A + 0(C), A| ei=o = A0 (1) + O(C 2 ), 
3 i^ =0 = 2^^Ti 1 )+0(C 2 ), (4.1) 

a 

the super-Hamiltonian constraint given in Eq.(2.6) gives rise to the following linearized 
elliptic equation for 

A0 (1) (r,oO - - 2 -^M^ V \r,a) + \ J2dlT^(r,a) + 0(C 2 ). (4.2) 

a 

2. The Equation for the Lapse Function 1 + 

In the family of 3-orthogonal gauges Eqs.(3.7) imply 3 K = 3 K^ = j^0(() ~ with 
F(i) (r, a) arbitrary numerical function of the same order. 
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To find the linearization of Eq.(2.12) for the lapse function we need the following result 



= -2 £ S 3 (a,ff % (r))m ( i V 7 = )(r,g) - 

w ) 2 + ^ E ^ - E ( 2 < - E ^ ^ 9 * c ^ ^w)) 

a ab a k^j 

m Qj da c(a, rfjir))] (r, a) - 2 £ ftft = 

j a 

= -2E T « " 2 E ^(2)' ^ rom ^9-(3-12), (4.3) 



where the following approximation was used 



'T,f)i(T)) 



^ro? c* + Ea q-2 ^ ia(r) - f A ±a ) 



2 



1 + 



+ / TrM ) AOO) + tt; o(0- (4-4) 

As a consequence, by using Eqs. (4.3) and (3.12), Eq.(2.12) becomes the following lin- 
earized elliptic equation for the lapse function 



An (1) (r,a) = -^ (1) (r,oO + i]Td a 2 ^ 

a a 

^ -ft •*„,(,-,?) + ^ + £ ig)(r,9). (4.5) 



5. TTte Solutions for <f>^ and 

The solutions, vanishing at spatial infinity, of Eqs. (4.2) and (4.5) for 0^) and are 

(i) 



{MfJ ] is given in Eq.(3.9)) 
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0(1) ( r >^) = 



c 



c 3 A 



+ E C (^rj-^ic^^r)))' 



m 2 c 2 



1^- 



i 

^/^([E(K) 2 -H-E^.^™) 
Eft^)^El])M)]- 



16n 



d 3 cti 



cr — 0"i 



<9 r g = 9 r ( i) +0(C 2 ), 



W-M = -e^[l + 2(r( 1 ) + 20 (1) )(r,a)], 



(4.6) 



n ( i)(r, a) 



r47rG 1 



,(17V) 



A 



G , yX 2 c 2 + E c («icW-%^c(r,«(r)))' 

1 + ^ L ; 

™ 2 c 2 + E c (^W-^icMiW))' 



+ / dVi 



3 d T 3 K(t, S x ) 



Ait \a — (Ti| 
^(i)rr(r,CT) = e[l + 2n ( i)(r,ff)]. 



(4.7) 



While 0(i) depends upon the tidal variables, the lapse function of these 3-orthogonal 
gauges is independent from them. However, while the volume 3-element 0m of the 3-space 
is independent from the inertial gauge variable 3 -^(i), the lapse function, connecting nearby 
instantaneous 3-space with different local York times, depends upon <9 r 3 /C(i) = d T j^ 3 K(iy 

At this order the spatial 4-metric 4 g(i) rs depends upon 2 <f>^ + : the first term describes 
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the instantaneous inertial part of the gravitational field in the 3-orthogonal gauges, while 
the tidal term the retarded one. Instead the component A gri) TT , relevant for local proper 
time, depends only upon the instantaneous inertial part of the gravitational field and upon 
the inertial gauge variable <9 T 3 /C(i). 

B. The Super-Momentum Constraints and the Shift Functions 

Let us now consider the determination of <7(i)( a )(b) |o^&> by means of the super-momentum 
constraints (2.7) and the determination of the shift functions from Eqs. (2.11). As said after 
Eq.(3.7), such a solution a (l)(a){b) \ aj A b implies vrf ) = ^ Y. a ^b e *ab (^a ) - T [ b 1} ) (r (l)(a){b) w 
in the York canonical basis of the linearized theory. 

1. The Equations for 0"(i)( a )(b) | a ^b 

Eqs. (2.7) for <7(i)( a )(b) |a^b takes the following linearized form {M^J is given in Eq.(3.10)) 

9 b a m(a)(b) (T, 3) ~\da 3 K (1) (r, B) + (Mf^J + 7aa 9 a H-a) {j, (?) = 

bj^a a 

= (d T d a ri 1} + ^d a {J2 d c n ( i)( c) ) - d 2 a n (1)(a) + 

c 

(4.8) 

where we have used Eq.(3.6) for U a and •faalab = o~ab — § (see before Eq.(2.3)). 

2. The Equations for the Shift Functions ^(i)( r ) 

Eqs. (2. 11) with a ^ b for the shift functions gives rise to the following linearized equations 

(d b n {1){a) (T,a) + 9 a n (1)(b) (r,CT))| a#b = 2o- (1)(o)(6) | 0#b (r,CT), a ^ b. (4.9) 

3. The Solutions for n(i)( a ) and °"(i)(a)(b) U^b 

By applying the operator <9& to Eqs. (4.9) and by summing over b, we get 

XI +<9a(9 6 n (1)(6) )](r,a)=2 X 9 6 (T(l)( a )(6)(r,?). (4.10) 

b^a b^a 

By putting Eqs. (4.8) for (7(i)( )(b)|a^b into Eq.(4.10), we get an equation containing only 
the shift functions 
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-[gCa *(i) + — ^— M (1)a -2C a n (1)(a) + 



(4.11) 



Eqs.(4.11) can be rewritten in the form (no more containing the condition b ^ a) 



&n {1){a) + -d a Q2 d b n {im ) {r,a) = 



(r,a). 



(4.12) 



If we apply J2 a d a to Eqs.(4.12) we get 



4 
3 

+ 



16nG 



(4.13) 



and this equation implies 



(4.14) 



As a consequence the final linearized equation (of elliptic type) for the shift functions is 



c 

+ ^(4^(9 a ri 1 ))-|^9 T (9 c 2 r«))](r ) a), (4.15) 



whose solution, vanishing at spatial infinity and depending upon the tidal variables, is 
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n(i)(a)(r, a) 



A Km + 



c 

4 G r«i«(r) - f A ±a (r, ^(r)) 



4ttG 1 



T,(7 



- E («ic(r)-^A ±c (r,7?;(r))) J 



d 3 ai 



\<7i - r}i{T) 



8lT |(7i — (T 2 | 



+ 



+ 



G 



r (fa i x - 

[4F as (r,^) (<(r^) - £ £ ^l^) 

n i * 

+ / ^ 



+ 



47T !<?! — <7 2 | 

F cs (r )( r 2 ) (vrKr,^) - £ 5 s " £ 
^(i)rr(T, a) = -en ( i )(r) (r,a). 



9c(o ; 2,r/i(r)) 



5 



a? /J' 



(4.16) 



Then the functions o"(i)( a )(6) \ a ^b have the the following expression 



V{l)(a){b)Wb 



\ (d a n (1)(6) +<9 6 n (1)(a) )(T,o ? ) = 



_o d a d b 
~ ~A~ 



8ttG 



3 *(D + -7s 



1 9„9 h 



a. 



J- t7a^fe kAUV) 

2 A ^ A • M W C 



+ 



+ s^(r« + r(')-i^|rw)= 
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= ~o 1^ ( 6 add b + 5 bd d a ) — > ^ ° — ( v, 

- J] («ic(r) - — A ±c (t,?7;(t))) / — — J-^ z —— 

V C /J All \a — (7i| |(Ti — 77i(T)| d 

/ r Q K-^(r))K-^(r)) ^ 



+ 



/ 47r |ff I 1 ffl | [ { 3R m + 2 rf) (r, ffx) + J Sa 2 



8ti Iff i — <7 2 



! + 



+ 



G 
2? 



d 3 ai 



\a — (Ji 



+ V d 16 d lc 

— ■*' I All |(7i — (7 2 

<i 3 <Ti 



«i 3 cr 2 



^ [4F h (T,?i) (ttKt,^) - S ° n E ^< 

n i 



<9c(ai,rfj(T))- 
dc(a 2 ,ffi(r)) 



+ 



+ 



9c(o ; 1 ,7/ i (r)) 



+ 



+ £ d 1(l «9 lc / 47r . ff 3 °! ff . ^(r.ffa) (<(r,«? 2 ) - £ 5- £ g i77 „ 



do? 
dc(a 2 ,ffi(r)) 



An <7i — <7 2 



(4.17) 



Both the functions ^(i)( a ) and <7(i)( a )(&) | a ^6 depend upon matter, upon the tidal variables 
and upon the spatial gradients of the inertial gauge function 3 /C(i) = ^ 3 -^(i)- 

Due to Eq.(4.16) we have that gravito-magnetism (described by 4 g^ rr ) in the 3- 
orthogonal gauges depends on both the instantaneous inertial and retarded parts of the 
gravitational field and upon the gauge variable <9 r 3 /C(i). 

The integral appearing in the shift function has the following expression 

f rf 3 ^ ( xac o K^WMH) ^ - 1 1 (sac K-<M)(^-^(t)A 

J 47r|o ! -5 ! i||ai-7? j (r)|3 \° ° |<?i-r?i(r)| 2 J " 2 \a-r}i(r)\ \° ' \S-m{r)\ 2 )' 

so that the contribution to gravito-magnetism coming from the mass current den- 
sity Mf$ r) {T,d) has the final form -™ £\ ^. o ( r ) _ & ^(r, j*(r)) + 



K-^(r)) 



«i(T)-^ i4 ± (T,*(T)) 



■(a-ifi(T)) 



\z-ni{T) 1 2 

From Eqs.(3.6) and (4.16) we get that the tidal momenta n s have the following expression 
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8ttG 



U- a (r, a) = d T R- a (r, a)-J2 Taa [d T ^ (4 lf> - 1 £ d\ lf>) + 

a c 

4tt G 1 , . „ 



c 

= ( E ^ s t -Rs - E s < 4 4 <«' - f S ft *<!5P) + 

6 a c 

M- a i = ^ - E ^ f ( 2 ^ - 5 E ^ f ) ' ( 4 - 18 ) 

a " "ft 

where we introduced the operator M^, which will be shown to be connected with the 
selection of the traceless-transverse (TT) part of the 3-metric on S r in Section VI. The 
relation between tidal momenta and tidal velocities depends on the inertial gauge variable 

C. The Contracted Bianchi Identities 

The contracted Bianchi identity (2.8) for d T <f>^ has the following linearized form 

d T (1) (r, a) = 6 d T q {1) (r, a) = l - ( E d a n {1){a) - 3 K (1) ) (r, a) + 0(( 2 ). (4.19) 



By using the solutions (4.6) and (4.16) for 0(i) and W(i)(o), this equation is identically 
satisfied at the lower order due to the conservation d^T^ = at the lowest order, see 
Eq.(3.13). 

The Bianchi identities (2.9) become 

d T n^0 + O(( 2 ), (4.20) 

consistently with Eq. (3.7) , which gives Trf\r, a) = + 0(( 2 ) in our 3-orthogonal gauges. 

It can also be checked with a lengthy calculation that also the contracted Bianchi identities 
(2.10) for (T(i)( a )(ft) \ a ^b are satisfied at this order. This check requires the use of Eqs. (4.8), 
(3.6), (4.9), (4.17), (4.6), (4.7), (3.13) (i.e. of the already found solutions of the linearized 
Hamilton equations) to get an expression which vanishes if we use Eqs. (6. 11), which are a 
byproduct of the linearized second order equations (6.4) for the tidal variables as shown in 
Section VI. 

As a consequence possible problems with the generalized gravitational Gribov problem, 
identified in Ref. [6], will appear at higher orders. 

Therefore the contracted Bianchi identities are identically satisfied at the lowest order. 
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D. The ADM Poincare' Generators 



The weak ADM Poincare' generators in the family of 3-orthogonal gauges are given in 
Eqs. (2.14) and (2.20), with the mass and momentum densities of Eqs.(2.4). They are the 
analogue of the internal Poincare' generators of Minkowski inertial rest-frame instant form 
[4]- 

By using Eq.(3.6) for the tidal momenta and Eqs. (2.4), (2.5) and (2.14), after a long but 
straightforward calculation (including various integrations by parts) we get the following 
form of the weak ADM energy at the second order 8 



l E ADM = Jd*a[M% V) + M% V) + 
c 3 



c 



+ IQ^G {^2^ R ~ a ~Yl 7aa d a ™(l)(a)) 2 + ^ °(l)(a){b) 
a a a^b 

- lCK (1) r-S {2) \ e , =0 \(r,a) + O(C 3 ) = 
= M (1) c+^Eadm(2) + McO(( 3 ), 

5(2) \e>=o = Yl [ 8 ( da ^(i)) 2 ~ 4 Yl 7sa da da Ra + 

a a 
+ J2 (2 Taa lla ~ fo) d a R a d a i%] + 0(( 3 ), 



M (1) c= J d 3 aM^ V \r, «?) = J> ^/ ™} c 2 + (/q(r) - ^ A ± (r, #(t))) * + 

^/^([E(w) : -(-«-E«-.^) 



+ 2c 



While the terms EJ a 7aa <9 a ^(i)( ) coming from the tidal momenta IT a are gravito '-magnetic potentials 
of the 3-orthogonal gauges (they depend on the shift functions), the term yj a ^ b cr ^i)( a )(t,) j determined 
by the super-momentum constraints, is a 3-coordinate dependent potential like <S( 2 ) (coming from the 
Gamma-Gamma term of the 3-curvature). 
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a 

+ ^ [i A4 (f7V) - A4 (f7V) - V A4 (C/y) -A4 (C/y) -- — 7U (c/y) Vl + 
c 3 U W A W L M (i)« A 4 A « a 7 J 

a a 

3 

+ 5 ^ ^ Rl + ^ ^ ^ 9a ) &) ' (4 ' 21) 

aft a 

Eqs. (4.6), (4.16) and (4.17) for (1) , fi(i)( a ) and <7(i)( )( 6 )| a ^ 6 have been used to get E ADM{2} . 
Since, as we will see in Section VII, the solution of the Hamilton equation for the tidal 
variables R a is proportional to ^, we see that all the terms not in the first line of E( 2 )adm 
are of order 4. 

Eq.(3.13) implies d T M (1) c = + ^ 0(( 2 ). 

In Eadm{2) we can see a kinetic term for the tidal variables R a , with the operator M n i 
appearing in Eq.(4.18) (the connected operator M ab appearing in the boosts (4.37) is defined 
in Eq.(6.4) of Section VI). Moreover, there is a potential for the tidal variables and bilinear 
terms in the matter. Finally, there is a term depending on the inertial gauge variable (non- 
local York time) 3 /C = ^ 3 -^(i) coupled to the divergence of the matter current density. 
However at this order the negative definite quadratic term in the local York time 3 Km = 
A 3 /C(i) appearing in Eq.(2.14) disappears due to the elimination of the tidal momenta with 
Eq.(4.18), the use of the solution of the constraints and of some of the Hamilton equations 
plus integrations by parts. As we will see in Eqs. (4.23) and (4.37) in the six Lorentz 
generators there is a dependence both on the local and non-local York time at the second 
order. 

For the other weak Poincare' generators the weak field approximation of Eqs. (2.20), by 
using Eqs. (4. 18) and by making integrations by parts, gives (all the terms in p r , 2 y and 

all the terms except the first in jJL are of order |§ on the solution for the tidal variables 
given in Section VII; the last term in the boost generators is a surface term which can be 
dropped with our boundary conditions) 

PADM = P r (l) +P r (2) + Mc0(( 3 )^0, 

P(U = / d 3 aM^ V r \r,a) = £ Vi (k^t) -|i ir (M,(r))) - 

i 

-- c Jd 3 aJ2 Frs(r, a) (<(r, B) - £ ^ £ Ql Vl ^fj^ ) , 
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f m = - J d 3 <j(^J2d r R- a M- ai d T R- b + 

ab 

a c 



(4.22) 



Jadm = J(i)+J( 2 )+ Mc °(C 3 )> 



3(1) 



= E ^ [vKt) (Kis(r) - 5l A ±s (t, a) ) - 17? (r) (« ir (r) - ^1 A ±r (r, a)) 

i 

- j k E (*i<*.*) - E «- E d -^0^) - 

u n i 

-a *ru{T,(T) \^±{r,(T) - > ^QiVi -q—^ 



3(2) 



J d 3 a (Mf x P (a r d s - a s d r ) 1 3 K (l) - 



a "' Vt ( 1 ) s A A x c 

<9 S . .(JJV\ 9 r 



+ 5 e <>:? ». - »* *) § 1 n» + 2 K>r | r?> - <>:> £ «•>) 

3 

■ g^[X> r 9.- ffS ^0 ^ M ^ 9 - ^ + 2 3 k {1) d r d s (r« - r«) 



+ 



+ 2 s r (r« + r?> - i E I ( r ?> " r ^)l ) ^ *), 

c 



(4.23) 



Jadm — —Jadm — 3(i) + 3(2) + McO(( 3 ) ps 0, 
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+ 



i/ A ,([E(kf-(-'-E«.^) 



+ 
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S) + 



^E^) + ^E(^) 

c a 

9 T r'"V 2 + ^F^fr« + r«"-^ 



ft 



+ 



+ 



J2 tfM 2 



c a 

— ^ [M ab d T T b ^ 

a,b c 



(4.24) 
(4.25) 

(4.26) 
(4.27) 
(4.28) 

A -(4.29) 



+ 



+ 



+ 



+ 



— ^ [M ab d T T b ^— ) + 

a,b c 
C 3 ^ dafrA / (1) (1) 1 ^ dlY ( P \ 9 a ft / s 4ttG ^ 

a^b c c 

c 3 ^r9 a 9 6 / 3 47rG^ ^jg; ^ 1 2 

n ^ d a M^ b dadbf AnG^ d c M^ 
2 2^ ~A A" I " ~ Z^^^J " 

c 

— L [(—5— J + — —J )( r '-) + 

c 



(4.30) 



+ 



E«) + (-) 

c 

(4.32) 
(4.33) 
(4.34) 
(4.35) 
(4.36) 

J L Z< Z_A 1UVI \ Z_A /J 

s c 

(4.37) 



We see that at the lowest order we get an unfaithful realization of the ten internal 
Poincare' generators M^c, p r ,^ ~ + 0(( 2 ), jTA, jT^ ~ + 0(£ 2 ) for the same matter 
in a rest-frame instant form approximately valid [modulo corrections at 0(( 2 )} in an ab- 
stract Minkowski space-time with the Wigner instantaneous 3-spaces coinciding with the 
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asymptotic inertial 3-spaces at spatial infinity of our space-times. They are the internal 
Poincare' generators of the system positive -energy charged scalar particles plus transverse 
electro-magnetic field in the radiation gauge studied in Refs.[12, 13]. 

Like in Minkowski space-time the internal 3-center of mass ff(r) of the isolated sys- 
tem 3-universe (i.e. "gravitational field plus particles plus electro-magnetic field" on 
the instantaneous 3-space S T ) and its conjugate momentum are eliminated by the con- 
straints p r {1) « + 0(C) and j^ d = - M (1)C 7f(r) w + 0(C)- They imply rf(r) = 
— M ^ - J d 3 a a r Ai^^\r, a) ps 0, i.e. the 3-center of mass of the 3-universe is in the origin 

of the 3-coordinates: the world-line of the time-like observer is the Fokker-Pryce 4-center of 
inertia F M (r) = 2^(r, 0) [4, 12, 13], to be interpreted as a version in the bulk of an asymptotic 
inertial observer. 

With the methods used in Refs.[4, 12, 13] it could be shown that in our asymptotically 
Minkowskian space-times there is a decoupled (Newton- Wigner) canonical 4-center of mass 
^com( r ) = ztl ( T jV( T )) °f the instantaneous 3-universes E T , carrying a pole-dipole structure 
(mass - Eadm and rest spin J^dm) and non-covariant with respect to an external asymptotic 
ADM Poincare' group with the same generators as in special relativity [4, 12, 13]. 



E. The Effective Hamiltonian in the 3-Orthogonal Gauge with a Given 3 K^(t, a) 

The restriction of the Hamilton equations of paper I to our family of 3-orthogonal gauges, 
given in Section II, can also be generated by using an effective Hamiltonian determined by 
the gauge-fixing procedure. 

Since our gauge fixings are 9 1 (t, a) ps and tt^t, a) — j^-^ F(t, a) ps with 3 K(r, a) pa 
F(t, a) an arbitrary numerical function, we must take into account the explicit r-dependence 
of this function which implies that the effective Hamiltonian is no more \ Eadm with the 
weak ADM energy given in Eq.(4.21). 

To find the effective Hamiltonian we should perform the following time-dependent canon- 
ical transformation: vr^(r, a) >->■ tt-(t, o) = 7r ^( r ! <?) — tJtg ^( r > ^ ( r > ^) = ^( r > w ^h 
all the other canonical variables fixed. Due to the explicit time-dependence the new effective 
Dirac Hamiltonian is H D = H D + ^ G J d 3 a 9F q T ^ (f)(T, a). In the new canonical basis the 
gauge fixings are 9 1 (t, a) ~ 0, vr^(r, a) ~ and have no explicit r-dependence. 

However to go to the reduced phase space we have to find the Dirac brackets and this 
would require the explicit solution of the super-Hamiltonian and super-momentum con- 
straints. This can be done in the linearized theory. 

From Eq.(4.18) we get (see Section VI for the inverse of the operator M n i) 



d T R- a = J2 M £ 



8ttG 



Tba 



4ttG 1 



01 



(4.38) 
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and this expression can be put in Eadm(2) of Eq. (4.21). As a consequence, the effective 
Hamiltonian of the linearized theory in the 3-orthogonal gauges is (use Eq.(4.6) for = 

1 + 60 (1) + O(C 2 )) 



Heff = M(i) C + - E A DM{2) + 



+ 



J d 3 ad T 3 K {1) (r,a) 



12nG 

= fd*o(M«P + M%P + 



6 

1 + A 



^^ + iE*r?>)(r,*)] = 



8ttG n { uv) ]_ m (uv) 



+ (E|^) S ^(i) + i^EE 9 ai2 B M a5 9 a i25 + 



+ — V r^n- + V^- -(45 M (uv) -^ T d M (UV) ) + ^ 3 K ( J) 



ab a 

-AnG 1 



d 2 



M. 



ab 

'8nG 



* 3A "<" I 1 + 1 ( " ^ M tP + i£ S » 2 r«") (r, *)] ) (r, *). 



a b 
,3 



+ 



127rG 



(4.39) 



As expected we get f (m (1) c 2 + £ iM (2)) = £ (m (1) c 2 + £ A dm( 2 )) + {M (1) c 2 + 

Eadm{2)i Heff} — (J^ acts only on 3 K^(t, a)). One could check that the Hamilton 
equations in the 3-orthogonal gauges are generated by this effective Hamiltonian 

In conclusion, when 3 K^(t, a) ~ F^(r,a) is r-dependent the effective Hamiltonian is 
not the energy and there are additional inertial effects like it happens in non-inertial frames 
in Minkowski space-time [4]. 
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V. THE LINEARIZED EQUATIONS FOR THE PARTICLES AND THE 
ELECTRO-MAGNETIC FIELD 

A. The Equations of Motion of the Particles 

Let us now study the weak field approximation of the Hamilton equations (2.18) for the 
particles. Since in this approximation we have 



1 »(i) ~ E c V^r) [n ( i)( c) + + 2 (1) ) (t)] 



l-€(r) 



:r,rf l (r)) + 0(C 2 ), 



(5.1) 



we get that the particle momenta (2.17) have the following expression in terms of the particle 
velocities 



Ki r { T ) 



rrii c 



Qi 



rrij 



A ±r (r,ffi(T)) + 



l-C(r) 



i)[(r) (l + 2(r( 1 ) + 2 (1) ) 



"(i) - E c i)f (r) [n (1)(c) + (lf> + 2 (1) ) #(t 



1-&) 



. + ^0(C) + ^0(c 2 ), 

^/ s rrii rrii 



l-^(r) 



or 



/t 7 > I T I 



TOj C 



TOj C 2 



Aj. r (r,77i(r)) + 



+ 



l-^(r) 



6 



-E*w(^i(^S?-|E*H l ? ) ) + 

c 6 

+ ^a| (4r "-5:f r ™ ) ) + 
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+ E«w>" ( - 4 -f z < V) + T '<> + 5 £ I r > ] ) + 

6 



+ i(9 r 3 ir (1) + ^§^9 r 3 ir (1) )l(r,rf l (r)) ) 
^ 1 -^(r) 



(5.2) 



where Eqs. (4.6), (4.7) and (4.16) were used to find the final expression. This result is 
implied by Eqs. (3. 11): the ultraviolet cutoff implies a small deviation from the expression 
in the free case. 

In the last expression the mass density the mass current density .M|^P and the 

tidal variables ri 1 ^, which will be shown to depend on the stress tensor TTA in Sections VI and 
VII, have to be evaluated by using the lowest order expression K ir {j) — ^ A± r (r, ^(t)) >->■ 

~2 



J^L (so that ^m^+(Hr)-^MTMr))) H- -7^=)- 



The particle momenta depend on the r- and spatial-derivatives of the inertial gauge 
variable 3 /C(i) = \ 3 -^(i)- If the York time 3 if(i) would depend also on the particle positions, 
we should make the following replacement d T 3 K^ h-> d T 3 K^\fj. + r)f(r) d b 3 K^. 

As a consequence, we can get the following second order form of the equations of motion 
(2.17) for the particles, implied by the Hamilton equations (we use the solutions (4.6), (4.7), 

(4.16) for (1) , n ( i), W(l)(r))> 



d_ 

^ dr 



\l 1 - Vi(T) 



(i + 2(r( 1 ) + 2 (1) )- 



HD " E c Vt(r) [n (1)(c) + (if } + 2 (1) ) r)f(r)] 



W(l)(r) 
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dr 



8ttG 1 



i + 2 r?)-^^™ + ££r?>- 



c 3 A 



A 



A 



C 3 A V— A 

c 



+ i§M4r«-£fr»W W (r«- ^ 



c 

^r+iEfOD 



+ 



+ 



= ^ == Ia (1) A ^ ( 1 ) r ~ A dc ^(De J + 



2 A 



to,- c 



2 F rs (r, ^(r)) ^ (r) - — — — + — — F ir (r, ^(r)), 



rriiC dr\\ 



TO,- C 



W[r,^(r)] = [d 3 aJ2 [(l + n (i) + 2 (r« - (1) )) W (n)a + 

+ (i-(r( 1 ) + 2 (1) ))n'„](T,a) = 

= / *° E [(i + *r z (^a'."' + E T w) + 2 rf' - 

a 6 

b 

+ (i + f^*™-i?>-i£§if>) w .]( Tl9 ), 



W (n)a (r,o ! ) = (2tt^(t,<7) - rijQjdacia^jiT))^ VkQkd a c(a,ff k (r)), 

j^i,k k^i 

W a (r,(?) = — E F ^( T '^) E VkQkd s c(a,fJ k (T)), 
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yi-^(r) a 

" -q^\(tMt)) = 

rriiC d . a /AnG 1 {U v) d a 



w [E #m i (^Sr 1 - 1 E ft + 



+ ^4 (4r « 1, -E| r i' ) ) + | 3 ^(.) + 

C 

. l!E± 1 M {uv) + r« + lv? r«l ^1 



(5.3) 



where W is the non-inertial Coulomb potential term and F ir are generalized inertial forces 
(now function of the inertial and tidal components of the gravitational field). As a conse- 
quence, the deviation from free motions is at the first order, consistently with the weak field 
approximation. 

Let us remark that in absence of the electro-magnetic field the final form of the equations 
of motion of particle % does not depend upon the mass Wj like it happens for test particles 
following geodesies. Therefore the masses are playing both the role of inertial and 
gravitational mass of the particles: their equality implies that the final form of the equations 
for particle i only depends on the masses rnj^i present in the mass density, in the mass current 
density and in the tidal variables. 

These equations of motion for dynamical (not test) scalar particles with a definite sign of 
the energy (implied by the Grassmann charges rji as shown in paper I) are not thought as the 
point limit of small extended objects as in Ref. [16] and do not contain terms corresponding 
to a gravitational self-force [16, 17] because rjf mf = (only terms rji r/j rrii rnj with i ^ j 
appear). As shown in Section VII with this description we can still get the energy balance 
when GW's are emitted. 

The Newtonian limit of the equations of motion for the particle in absence of the electro- 
magnetic field will be studied in the third paper [18]. 

B. The Equations of Motion of the Transverse Electro-Magnetic Field 

The weak field limit of Eqs.(2.19) is 
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d T A ±r (T, &) = S rn ^V) + "(1) " 2 - 0(1))] < + 

na 



+ n( 1 )-2(ri 1 ) + 2 0( 1) ) + 



y/m^ + E b far) - % A ±b )' 

E fe (ri 1) +20( 1) ) (K lb {r) - ^ A ±b )\ 
7 7, \2 — j ~ "WW 



b 

[d b 0(i) - d b n ( i) + <9 6 (r« + if } )] F a6 ) - ^ n ( i) (6) <9 6 < 



+ 2^ ^ E l^ n (i)W Q^-a d b n [1){a) 



0a b 



(5.4) 



The explicit expression of these equations in terms of the instantaneous inertial and re- 
tarded gravitational quantities can be obtained by using Eqs. (4.6), (4.7) and (4.16). At the 
lowest order and using Eq.(5.2) these equations imply the special relativistic wave equation 

DA ±r (r, a)= Y.nsm SrnP± s (°)<>s m J2i Vi Qi vT( T ) <* 3 (^> Vi( T )) , whose Lienard-Wiechert so- 
lution was found and put in Hamiltonian form in Ref. [15]. 

The first of Eqs. (5.4) can be inverted by iteration and the resulting form of the transverse 
electro-magnetic momenta is 



*1(t, 3) = E P T^) ([1 - n ( i) + 2 (r« - 0(i))] d T A ±a - 

a 



(5.5) 



If we put Eq.(5.5) into the second of Eqs. (5.4) we get a wave equation for Aj_ r (r, a) con- 
taining extra terms in d T A± a and d b d T A± a , which will be studied elsewhere. By using the 
solutions 0(i), ri(i), n(i)(a), we see that this equation has as sources both the matter and the 
tidal variables. Since in the next two Sections we will show that the tidal variables are 
determined by the stress tensor TTA in a retarded way, the final equations for the transverse 
electro-magnetic field depend on the matter both in an (action-at-a-distance) instantaneous 
way and in a retarded way. If one could find a Lienard-Wiechert-type solution of these equa- 
tions, the final form of the particle equations (5.3) would be of coupled integro-differential 
equations instead of differential equations like it happens in absence of the gravitational 
field. 



53 



VI. THE LINEARIZED SECOND ORDER EQUATIONS FOR THE TIDAL VARI- 
ABLES R a IN THE 3-ORTHOGONAL GAUGES 



In Section IV we solved the equations of elliptic type for 4>{i), n^, 7i(i)( r ) and <7(i)( a )(b) \ a ^b- 
The solutions depend on the tidal variables R a . 

We must now study the linearization of the second order equations (2.16) for the tidal 
variables by using the solutions of Section IV. We will see that also in these non-harmonic 
3-orthogonal gauges we get wave equations, but they will also contain the information that 
the final 3-metric is traceless and transverse (TT). 

A. The Linearization of Eqs. (2. 16). 

For the three integrals appearing in Eqs. (2. 16) we get the following linearization: 
a) the last integral in Eq.(2.16)) becomes 



5M(t, (j\ 



/«M 1 + n < T '*V 5R- a (r,v) 

, J2a laa Uia(r) - ^ A ±a j 

= -J2^m(r))m( i V J = )(r,g) + 

+1 [ Yl 7«« W) 2 - E Too F 2 ah - J2 Taa (2<-E Vk Qk d a c(a, %(r))) 



c 

a ab a kjtj 

-2\ 



VjQjdac(a,rf j (T))^(T,a) + 0(mc( 2 

3 

= -Y^laaT^ + 0(( 2 ), from Eg. (3. 12). (6.1) 

a 

b) the first integral in Eq.(2.16) becomes 



/ 



dTOi [1 + n(r, <Ti) 



=o — 



<Ji2a(r, <t) 

= ^T, ^ 0(D (r, ?)-2^(2 7 aa 75a - hi) $1 R~ b (r, B) + 0(( 2 ). (6.2) 



ab 



c) the second integral in Eq.(2.16) becomes 

'^^aO ^^ l^o = 2^ T ^n (1) (r lff l+0(( 2 ), (6.3) 
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As a consequence, the linearization of the second order equations (2.16) for the tidal variables 
Ra is 



d 2 T R- a (r,a) = AR- a (r,a) + Y,Ta 



d T d a ra ( i)( a ) + 



+ dln {l) +2dl<P w -2dlY^ + 



8nG 



^ by using Eqs.(4.6), (4.7), (4.16) 



nY,M- ai R- b = E- a , 

1 

M s5 = % - 7 - a V 2 75a ~ 2 ^ 756 A") ' 

a 6 

4 ^E^4^(4< ) r ) -|E^< ) r ) ) 



+ 



<9 : 

I r) 71(10 1 a \ r rpbb 

(1) A ^ (1) 

b 



4 



b a 

M ab = TaaTbbM- ai = <U (l - 2 I) + I (l + I) |, 



^ M ab = 0, il4 5 = ^ laa Tbb M ab , 



(6.4) 



ab 



where we used Eqs.(4.6), (4.7), (4.16) and ]T S Taalab = <U - §, E a 7aa75a = E a 7aa = 
0. 

Therefore we get the massless wave equation (with the flat d'Alambertian □ associated 
to the asymptotic Minkowski metric) not for the tidal variable R a but for the quantity 
J2b M a iRi, where M n i is the spatial operator already found in Eq.(4.18). 

B. The Meaning of the Operators M- ab and M ab 

Let us show that the operators M ab - and M ab are present to select the TT part of the 
spatial metric 4 g(i) rs = — e 3 g(i) rs = — e^rs + 4 ^(i)rs = — e^rs + 2(1^ + 2 0(x))j, namely 
that they define the polarization pattern of the gravitational waves in these non-harmonic 
3-orthogonal gauges. 
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In Ref.[20] it is shown that in every gauge we can make the following decomposition of 

h(i) rs (r, a) 



4 h {1)rs = A hf^ rs + -5 rs H w + ^ (d r e(i) s + d s €(i) r ) + (d r d s - ^5 rs A) (6.5) 



1 
2 



1 

3 



with Y^ r ®r £(i)r = and 4 hJ^ rs traceless and transverse, i.e. J2 r 4 ^(i>r = °> J2 r ®r A hf^ rs 
0. The functions #(i), A(i) and e^iy have the following expression 



H {1) 
A(i) 

C(l) r 



3 1 

2 A2 



$3 (ft, 9, - 1 <U A) % )fl0 = -3 6 £ || it 1 ), 



<9 r 



^ ^ v u n {1)ur 



= -^i(r? , -E|r?>). 



4/, _ 4tTT 

ft(l)rs - ft(l) rs _ 6 



(^u>-£§<K,+ 
+ 2^(r« + r«-i^|r«). 



(6.6) 



In the last line of this equation we have given the form of the 3-metric implied by Eq.(6.5). 

As a consequence the TT part of the spatial metric is independent from 4>{i) and has the 
expression 



AhT (l)rs = 4/i (l)rs + -Z Srs ( ^3 ^ 



A 



(l)uv 



i a. 9., 



+ 2 ^ £ ( 53 Ah ^ uv + 53 4/l ( 1 )««) ~ 53 X ( 5r 4/l W us + ^ 4 ^( 1 ) ur ) = 



<9 U 



= — e 



[(2 r« + £ I r«) - 2 ^ (if> + r?>) + ^ £ § ^ 

u u 
^ ^ 'Prsuv h(l)uvi 



d r d x 



6„ 



A J 2 



\ ( 6rs + 



d r d s \ d u d v 



A J A 



(5 rv d s + 5 SV d r ) + (5 ru d s + 5 SU d r ) 



(6.7) 
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where V TSUV is the projector extracting the TT part of the spatial metric. It satis- 
fies y~^ fT , Prsuv Puvmn Prsmn Psrmn Prsnmi &r Prsuv ^ u Prsuv Oj 
Prruv Prsuu 0. 

With plane waves, i.e. 4 /i(i) rs (r, a) = A TS {r) e in ' a + cc, we recover the standard projector 
Vrsuv ^rsuv = P r u P S v — \ Prs P U v with P rs = 5 rs — n r n s (n is a unit vector orthogonal to 
the wave-front). 

The diagonal elements 4 hJ^ aa of the TT 3-metric contain the operator M a b of Eq.(6.4) 
since we have 



X^ = -2 e [(i-2|)r« + i(i + |)i:|r« 

c 

= -2e^M ac r«, 

c 

=► □ 4 / i fl)aa=-26^7,a^. (6.8) 

a 

If we apply the decomposition (6.5) to the stress tensor TjT^ of Eqs.(3.12), we get 
1 1 1 

IJl) = g #(1) ^rs + T(i) rs + ^ ^(i) s + ^ s + (^r ^ s _ 3 ^) 'V 1 )' 

E^-E^r'-E^-". 

r r r 

r 

rs 

e (1)r = 2^(9 u T ( ^-i9 r T ( «p-2|^(9 u ^-i^A)T ( ^ 

M UV 

rp{TT)ab _ \ ^ rped 
J (l) - / , ^o6crf-t(i), 

T (i) = 2^ [ 5ac 5ad ~ 2 "X" 2 A" ~A~ ~ 2 cd ^ ~ A" J 

(6.9) 

where T^ T ^ ab is the TT part of the stress tensor. 



By using Eqs.(3.13), i.e. «9 T ^^ y) = -£ c ^^+3^^, d T .Mf y) =-£ c <9 C T ( - + 

<9 A K$, so that = Y. cd d c d d Tf* + d A (d T Tlf^ + £ c d c K$), we get the following 

expression for the source term appearing in the second member of Eq.(6.4) 
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V -v- E- — — T {TT)aa E- - — V t- T {TT)aa 

7aa £/ ° ~~ r 3 1 (1) ' - C/a ~ r 3 > aa J (1) 



(6.10) 



Therefore we get that the TT part of the 3-metric satisfies the wave equation 



= -2 e $> a ^ = 



167TG rp(TT)aa 
r 3 1 (1) 



4 



n 4 /)-- - — t 



167T G rp(TT)ab 



(1) 



(6.11) 



The second line of Eqs.(6.11) derives from the first line due to the transversality property 



TT 
(l)aa 



+ 



of %? rs and T$;i which implies 9 b (p^ ba + e^T^) = -d a {u±K 

c WnG T (TT)\ o_ 
6 C 3 I (l)aa J - U - 

Therefore, even if we are not in a harmonic gauge, the Hamilton equations in the 3- 
orthogonal gauges imply the massless wave equation for the TT part of the spatial metric. 

Once Eq.(6.11) is solved for 4 /i[^ = -2e M ab Y { ^ \ see Eq.(6.8), in terms of T^ T)aa , 

we can find the solution for R a = XL laa if we succeed to invert the operator M a j,. 
Then, if we put this solution in the last line of Eqs. (6.6), we can find the spatial metric 
4 /i(i) rs in the 3-orthogonal gauges. 



C. A Generalized TT Gauge 

At this stage the 4-metric of Eqs. (3.2) and (3.5) has the following form 



/ 2n (1 



9(i)ab = Vab + e 
= 4 Vab + e 

+ o(c 2 ), 



(i) 



-n ( i) (r ) 



\ 



^ -2 (r T (1) + 2 (1) j 5 rs y 



+ o(C 2 ) = 



-f 3 K (1) + A r (ri 1} ) + ^{matter) \ 



\-^ 3 K {1) +A s (T ( a L> ) + P s (matter) B r (T { a>) + >y (matter) 5 rs J 



+ 
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A-(rf>) = 
s r (rW) = 

a(matter) = 
6 r (matter) = 
j (matter) = 



-^IK'-Ef^). 

c 

-2(r?> + iEf^»). 

c 



47rG 1 



_ d M (uv)\ 



AM) 



8ttG 1 (c/y) 



(6.12) 



by using the expression of <f>m, nm and nmM given in Eqs. (4.6), (4.7) and (4.16), re- 
spectively. We have explicitly shown the dependence upon the inertial gauge variable 
3 /C(i) = ^ 3 i^(i), the instantaneous inertial effects a(matter), (3 r (matter), ^(matter) and 



the retarded tidal effects ^(ri 1 ^), B r (T 

Let us consider the following coordinate transformation on the 3-space S T (endorsing it 
with r-dependent new radar 3-coordinates) 



f = r , a r = a r -^ r {1) (T,a), 

1 d r /,„ m ^ dl 



*(D M = -ig (4 I?> -E| r«) (r, a) = 0(C), 



(6.13) 



and the associated new 4- metric 4 g(i)AB( r ) <?) = ffx ffs- 4 fi , (i)cz?(' r ) o 2 ) replacing the one of 
Eq.(6.12). The inverse transformation is r = f , a r = a r + x I / (i)(r, <r) + 0(£ 2 )- 



Since we have = «-f^, ^ = -^5, we get ffl = S^^gfl+O^), 

^pf = — {1 Q T T,a + 0(( 2 ). In particular, by using Eqs.(4.6), (4.7) and (4.16), we get (since 
f = r and <t, we always use r) 



i 9(i)T T (r,o-) 



= e 



g (l)TT (r, a) + 0(( 2 ) = e [l + 2 n (1) ] (r, + 0(( 2 ) 
2 ^ 3 A' (1 , + 



4- 



0(i)tt(t, cr) 



c 3 A 
9* ( i)(r, a) 



9r 



+ o(c 2 ) 



<9^ 



(n ( i )W + ^)M) + 0(C 2 ) = 



^ 3 

A 



9r 
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%i)rs(r, a) = 4 g rs (r, a) - e (d r ttfo + d s ttfo) (r, B) + 0(( 2 



= -e[S r 



(5 rs [l + 2 (r« + 2 (r, <?) + [a r ttfo + & *(i)] ( r > ^)) + °(C 2 ) = 



= Vrs + 



4 7 TT 



(l)rs 



87r G o~ rs >(c/y)' 
c 3 A^W . 



M) + 0(C 2 ). 



(6.14) 



In the last line we used Eq.(4.6) for 0(d to recover the TT 3-metric hT^ rs of Eq.(6.7). As 
a consequence it turns out that the new 4-metric in the new radar 4-coordinates has the 
following expression 



4- 4 i 

9ab = Vab + e 



/ -2 ^ 3 fsT ( D + a (matter) -J 3 fT ( D + ^(matter) \ 
y 3 K(d + (3 s (matter) e 4 hj^ rs + <5 rs ^(matter) j 



o{C). 



(6.15) 



Therefore the coordinate transformation (6.13) leads to a generalized TT-gauge whose 3- 
metric is not 3-orthogonal due to the presence of the TT 3-metric. Also in absence of matter 
it differs from the usual harmonic ones, whose instantaneous 3-spaces are all Euclidean, for 
the non-spatial terms depending upon the inertial gauge variable 3 /C(i) = \ 3 K(i)(r, a) (the 
HPM form of the gauge freedom in clock synchronization). 

If the matter sources have a compact support and if the matter terms \J^fi^ and 
~k "Mrfy? are negligible in the radiation zone far away from the sources, then Eq.(6.15) gives 
a spatial TT-gauge with still the explicit dependence on the inertial gauge variable 3 /C(d 
(non existing in Newtonian gravity), which together with matter and the tidal variables, 
determines the non-Euclidean nature of the instantaneous 3-spaces. 



D. Inversion of the Operator M ab and the Resulting Form of the Second Order 
Equations for R a . 

Since Eq.(6.8) gives J2 b M ab Y [ b 1] = -f A hJ^ )aa with the diagonal elements of the TT 3- 

metric satisfying the wave equation (6.11), i.e. D 4 hJ-^ aa = — e T^ T - )aa , to find the tidal 

variables R a = ^2 a 7s a associated to one solution of the wave equation we have to invert 
the operator M ab . 

If we introduce the functions 



_ t4,TT 
9a — — 2 n (l)aa 



E = ^ + -A 1 + f ) E - 2 x r i : 



A / ^— ' A 

c 



E n i-i o 87T G (TT)aa 

g a = 0, ng a = —£-Tt 1) , 



(6.16) 
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we can write 



W =9a-da*h-\Y l d ° i *h> 



where is the function generating the coordinate transformation (6.13). 
If we define the quantities 



_ ivTT I 
Z a b\a^b ~ n (l)ab\a=£b 



d a d b 



(r?> + r?>-±X;§r?>) 



2(d a *\ 1) + d b W{ 1) ^ b , (a^b), 



Z a b\a^b — Z ab [gc\\a^b- 



(6.17) 



(6.18) 



they turn out to be functionals of the g^s expressible in terms of the functions of 
Eq.(6.13). Therefore, also the functions can be viewed as functionals VP®^ [<7 r ] of the 

9aS. 

As a consequence, the TT 3-metric of Eq.(6.7) can be written in the form 



A h TT - 9 f 
n (l)ab — ~ At 



2A ^ A 



5i ^12 [g] z 13 [g] 

Zl 2 [g] 52 ^23 [5] 
*L3 M ^23 M 53 



(6.19) 



The transversality property of the TT metric, shown in Eq.(6.18), allows to express z a b[g c 
as the following functional of g a 



^23 [5c] ( T > &) 

^13 [5c] (t, B) 

zu[g c }(T, B) 



2 3 

\ f ^ da 2 da 3 [-dlg 3 -d 2 2 g 2 + dl gi } (r,a\a 2 ,a% 

1 3 

\ f ^ da 1 da 3 [-dl 9l - 81 g 3 + d\ g 2 ] (r, a\ a 2 , a 3 ), 

1 2 

\ f da 1 da 2 [-d 2 2 g 2 -d 2 1 g 1 + d 2 g 3 \ (r, a 1 , a 2 , a 3 ). 



(6.20) 
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By using J2 a 9a = 0> Eq(6.20) can be put in the following form 

z ab [g c }\aAT,t) = -l da a da b [(A - d 2 a ) g b + (A- dl) g a ] (r , a a , a b , a c ^ a,b ). (6.21) 

Z Ja^b 



Due to Eqs.(6.18), the functionals [g c ] of the <? c 's are related to the functional given 
in Eqs.(6.21) by the equations 

(d a * b {1) \g c ] +d b y a {1) [g c ]}\^ b = z ab [g c ]\ a ^ b . (6.22) 

The solution [g c ] of this equation is (e 2 ,^ is symmetric in e and / and vanishes if a = e 
or a = f) 

1 f a2,ai ( \ 

^li)[gc](r, a) = g / da 2 da 3 ^d 3 z 12 [g c ] + d 2 z 13 [g c ]-d 1 z 23 [g c ]J(T,a 1 ,a 2 ,a 3 ), 

1 3 

*(i)[^c](t, <?) = - / do- 1 do- 3 ^ 3 2 23 [# c ] + d 2 z 12 [g c ] - d 1 z 12 [g c ]J (r, or 1 , a 2 , o- 3 ), 
^?i)[^ c ]( r > d) = \ J da 1 do 2 [ d i z 23[g c ] + d 3 z 12 [g c ] - d 2 z 13 [g c ]^j (r, a 1 , a 2 , a 3 ), 



or 



T e .rrf 



^(i)[^](r,a) = - ^{taeff J rfa 6 ^ (^<9 e 2 /a [# c ] + d f z ea [g c ] - d a z ef [g c ]j (r, a a , a e , a f ). 

ef 

(6.23) 

By using Eq.(6.21) and J2 a ga = we get the following form of Eqs.(6.23) 



( f \2 rcr^ ,<JJ r&,a' per" 

*(l)[ffc](T,5) = - Y.-^f- / / d&Cd&f / ^ 

e f J J J 

[(A - (9 2 )A g f + (A- d 2 )A g e - d 2 e d) (g e + g f )] (r, ~a\ o- e , a*). (6.24) 



Then by using Eqs.(6.17) we get 

r« = g a -d a ttfo [g b ] -lY, 9 * *<U W = 

c 

= E Kb M bc r?) = £ M", 1 ff6 , (6.25) 



6c 
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and this is a definition of the inverse operator M ah x by means of Eq.(6.24). 

If g c = p c is a solution of the wave equation □ g c = T^ T ^ CC (see Eqs.(6.11) and (6.16)), 
then the associated tidal variables (our physical radiative degrees of freedom replacing the 
TT 3-metric of the standard approach in harmonic gauges) are 

R-a = Yl 7*. T a ] = E 7a- M£ Pb = J2 7*. {pa - d a [ Pb ] - \ £ d c ¥fo [p b ]) , (6.26) 

a ab a c 

with the functional ^ a {1) [p b ] of Eqs.(6.23). Eqs.(6.4), (6.10), (6.11) imply that these tidal 
variables are solutions of the wave equations 



(1) o_ 8n G (TT)aa 



(TT)aa 
) 



(6.27) 



In the next Section we will study the solutions for the tidal variables, because they are the 
HPM gravitational waves with asymptotic background in our family of 3-orthogonal gauges 
with non-Euclidean 3-spaces E T . 

From Eqs.(2.3), (3.6), and by using ^aal&b = <5ab — §, we get that the extrinsic 
curvature tensor of our 3-spaces in our family of 3-orthogonal gauges is the following first 
order quantity 



3 K(i) rs = <7(i)( r )( s )| r ^ s + S rs y- 3 K(i) - drT^ +d r ri(i)(r) — 9 ri(i)( )j, 

a 

(6.28) 

with fi(\)( r ) and cr(i)( r )(s)\r^s given in Eqs.(4.16) and (4.17), respectively (they depend on 
3 /C(i) = j^ 3 K). Therefore, our (dynamically determined) 3-spaces have a first order devia- 
tion from Euclidean 3-spaces, embedded in the asymptotically flat space-time, determined 
by both instantaneous inertial matter effects and retarded tidal ones. Moreover the inertial 
gauge variable 3 if(i) (non existing in Newtonian gravity) is still free. From Eqs.(2.5) we see 

that the intrinsic 3-curvature of these non-Euclidean 3-spaces is z R\e* = o = 2 ^ a d^T^: it 
is determined only by the tidal variables, i.e. by the HPM gravitational waves propagating 
inside these 3-spaces. 

If we use the coordinate system of Eqs. (6.13) (6.15) to go in the generalized TT gauge, 
we can introduce the standard polarization pattern of gravitational waves for A h^J (see 
Refs. [20-22]) and then the inverse transformation allows to rewrite it in our family of 3- 
orthogonal gauges. 
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VII. POST-MINKOWSKIAN GRAVITATIONAL WAVES WITH ASYMPTOTIC 
BACKGROUND IN THE 3-ORTHOGONAL GAUGES 

In this Section we study the solutions of the linearized equations for the tidal variables 
Ra(r, a) and the TT 3-metric 4 hJ^ rs (r, <?), namely the PM gravitational waves with asymp- 
totic background in the family of 3-orthogonal gauges. When needed we assume the validity 
of multipolar expansion of the energy-momentum tensor, which is reviewed in Appendix B. 



A. The Retarded Solution 

Since in Eqs.(6.27) we have the flat wave operator □ = d 2 T — A associated with the 
asymptotic Minkowski 4-metric, we give the solution as a retarded integral over the past flat 
null cone attached to the point (r, a) on the instantaneous 3-space S T at time r by using 

the retarded Green function G(t, <t;t' ',<?') = — 6{t — t) S ^ T ~^f_^ - (□ G(r, a; r', a) = 

8(t'-t) 5 3 (a-a)). 

The retarded solution of Eqs.(6.27) is 

R- a (r, a) = J2 Taa r?>(r, <?) = £ 7aa M~\r, a) (F b T ™(r, " 

a ab 

2G ^gVM^^ 

J CT- (7i / 



C 3 



r« = ^ 7saJ R s , (7.i) 

a 

where Fi feom ' ) (r, c?) is a homogeneous solution of the flat wave operator (DFa hom ^ = 0) and 
the non-local operator is defined by means of Eq.(6.26). 

The condition of no-incoming radiation is Fa hom \T,a) = 0: it uses the flat light-cone of 
the asymptotic Minkowski metric at r — > — oo. Therefore there are only outgoing gravita- 
tional waves. 

With our matter (point particles plus the electro-magnetic field) we do not need to 
solve the Hamilton equations independently outside and inside the matter sources with the 
subsequent matching of the solutions like it is usually done with compact matter sources. 
Even if j^ T ) aa j s assumed to have compact support, what is relevant in Eqs.(7.1) is the 

behavior far from the support of the non-local quantity determined by applied to the 
retarded integral. 



Since we have T^ ab {r,a) = V abcd T^(T,a), with the operator V a bcd = Vabcd(&) defined 
in Eq.(6.7), the chosen solution (7.1) has the following form 
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6(t-\ 


a — <J\\ 


-n]) 




\a — o\ 





2G ~ f 

R- a (r,a) = -—Y^laaM- b l {r,a) J dn d 3 a, 9 (r - n ) 

ab 

Vbbuv((?i) T^(n, cti) = 

C° * ' J (7 — (7l \* ' W / T-tT-\<T-(Tl\ 



ab uv 



Analogously the solution for the TT 3-metric ti^ ab , satisfying the wave equation (6.11), 



is 



\ C r \ Y^uv 'Prsuv(^l) T^(t, <Ti)\ 

h (i)rs\ T i a )= e 3" d a l — — . (7.3) 

c J I a — G\ I 

Let us remark that, since we have DV a b c d = Vabcd^, another solution of Eqs.(6.11) is 
(in general it is a solution differing from Eq.(7.3) by a homogeneous solution of the wave 
equation) 

'hV»(r, S) = -e '4 £ P_(5) / ^ HCMfzgLkg) . (7 . 4) 

v y C J z — ' / (7 — <7i 

Ml) 11 

To compare Eq.(7.4) with Eq.(7.3) we need the following integral representation of the 
operator V rsuv {a) 9 

fl T (r,a) = J2 V rsu^)f UV (r,ap J d 3 a, ^^-^•/"(^l), 

uv ' uv 

d>rluv(<? - <?l) = ^ ((^ ru ^ + ^™ ^« ~~ ^ rs ^ uv ) ~~ ^) + 

ra $sb <^m>J 



+ 

ab 



5 ab 


a — U\\ 


2 _ ■ 


3 (a a - 


af) (a b - a\) 


47T 


a — u\\ 


5 



+ 



+ 5 rs / d 3 a 2 



5 rs 


a - a 2 \ 


2 _ • 


3 (a r - 




47r \a — a 2 \ 


5 



guv 


o 2 — 3\, 


2 _ 


3 K - 




4:71 \a 2 — <?i\ 


5 



)■ 

(7.5) 



Wc have = -5 3 (a), d r \5 - 5i\ = d r = d r d sW ^ 

" T^r [) . dr f(r -\a-a 1 \) = -^d T f(r - |ff - ad)- 
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where the integral kernel d^ cd is a function only of the difference a — B\. 

As a consequence, the two retarded solutions (7.3) and (7.4) for the TT 3-metric with 
the no-incoming radiation condition take the form 



I .. I 



A7TT , - 4G /■ 3 V-^T,^ -» \ /" j3 T (l)(^~ I <?2 ~ Ol |,Ol) 

tyi) <t) = -e — / <To 2 2^ d rsra (a - ^2) / dVi — — — , (7.7) 

c d y ^— ' J 0-2-0-1 

respectively. To make the comparison either an explicit form of the energy-momentum tensor 
or its multipolar expansion is needed. 



B. PM Gravitational Waves with Asymptotic Background: Multipolar Expansion 
and the Quadrupole Emission Formula 

To look for the PM relativistic quadrupole formula we need the multipolar expansion of 
the energy-momentum tensor in terms of relativistic Dixon multipoles expressed in our rest- 
frame instant form of dynamics. In Appendix B there is a review of such multipoles based 
upon Ref.[23]. See chapter 3 of Ref. [21] for a review of the standard multipole expansions 
used in the literature (see Refs. [24, 25]). 

1. The multipolar expansion for the tidal variables R a 

By using the multipolar expansion (B2) centered on the center of energy w^(r) = 2 M (r, 0) 
(r/(r) = 0) and the operator V rsuv defined in Eq.(6.7), from Eq.(7.2) we get by making 
integrations by parts (assumed valid) 



/\ 3 T<j^»(r-|*-*i|;*i) 
a o~i 



c 3 j \a — G\ 



c 3 nl 



n=0 uvr\..r n 



q ri..r n \m 


>(t-< 


0- Oi|) 




— Oi 





[o au o av 2 {L w+ 2 ( + Aj A, 

(S 9lu \r> \ ^"^(Ol) 
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2G 



CO / \ n 

E ^- E a---** 



c J 14 — ' ni 

n=0 n..r. 



qn..r n \aa^j. _ 






— * — * 





) k=0 



5E( 



9i ri ...9i rt] 



qn..r n \uu^ T _ 


<7-<7i|) 




—* —* 

a — <j\ 





) 1 5*1=0 



+ 



qn..r n \uvf T _ 


o-e\\) 




a — B\ 
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The n = term has the following expression 
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(7.9) 



In the last expression we used Eq.(B13) and (B14) and Y^ r V rrU v = 0. The last line used 
Eq.(6.26). 
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The first term of Eq.^7.9) gives the standard quadrupole emission formula modified by 
the non-local operator M~ b x of Eq.(6.26). The higher terms would give the contribution of 
the mass octupole, momentum quadrupole and so on (see Section 3.4 of Ref. [21]). 



2. The multipolar expansion for the TT 3-metric 4 /i(i )rs 

By putting in the solutions (7.6) and (7.7) the multipolar expansion (B2) with rf{r) = 
we get 

An 00 

4.TT 



qn...r n \uv^ T _ 


o 


-*i|) 




a — 0\\ 





n=0 r\..r n uv 

-dZM)- (7-10) 



da?...da\ 



<h% r .(T,#) = -e^± L2L [ d *o 2 J2 <fZ v {a-B 2 ) J2 

C n ' L - J „.„. „ 
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fqri-r n \uv^ T _ 




^|) 


do r 2 \..do r 2 n 1 


K 


°2 





n=0 " uv r\..r n 

= %] rs (r,o). (7 

These two expressions can be shown to coincide with the change of variable a — 3\ = o 2 
(<fVi = d 3 o 2 ) and by making suitable integrations by parts (we assume that the integrations 
by parts can be done). 

We shall use Eq.(7.11), for whose manipulation we need the formula 





fqn...r n \uv^ T _ 




**\) 


do r 2 \..do r 2 n 1 


K 


&2 





^ n\ d k ( 1 \ d n ~ k q ri - r ^ uv (T- \ o 2 |) 



k=0 



k\{n-k)\ do? ...da r 2 k \\a 2 \J da 2 k+1 ...da 



'2 



™ Tl ' f) k { 1 \ 

_ \ r ( \n-k _ I I an-k n...r„\uv t 

" t^ k\{n-k)\ K ' da r 2 \..da r 2 k \\ a 2 \) T q 1 



T~ I 0-2 \)n {2 )r k+1 -n {2)rn , 



O On 

n r = -r-TT, n (2)r = ——. (7.12) 
I o I I o 2 I 

Therefore Eq.(7.11) can be put in the form 

00 n l \2n-k 

Ar/TT 



n=0 fc=0 7 uv 

E a^W fe) arVl rn ' w(r - 1 " 2 



(7.13) 
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To study the behavior of Eq.(7.13) at big distances, i.e. r = \a\ » 1, we use the 
following results shown in Appendix B of Ref . [26] 

rv) * ^- =► f TT rs(ff) = / dv 2 ^ - - 2 ) ^ ^ + o{i/r% 

uv 

r V) ^ ^ / TT -(-) = / ^2 E - *o + o(i/r 3 ), 

Ml) 

rW H- — =► / TT rs (<?) = / dV 2 £ - 30 f V2) H- 0(l/r 3 ), n > 3. 

Ml) 

(7.14) 

As a consequence at great distance only the terms with k = give the dominant contri- 
bution in Eq.(7.13) 

A ti$ ra (T,3) = -6— I d 3 °2 d rL(B~ft) 

n=0 uv 

an ri ...r n \uv( I ^ |\ 

E »(2)n-n (2 )r„ ^ ,i| ' 2l) + 0(l/r 2 ). (7.15) 

n..r„ 1 1 

As shown in Refs. [20] (see also Ref. [25]) we have the following result 

r v) = r(\ b i), 

f T rs(ff) = [ d 3 a 2 CL(<? ~ ft) f UV (\ ft I) = E A rsuv{n) r(\ B I) + 0(l/r 2 ), 

(7.16) 

where A abcd (n) is the algebraic projector (defined after Eq.(6.7) for plane wave solutions) 

1 CX^ 

Kbcd{n) = {5 ac -n a n c )(5 bd -n b n c ) - - {S ab - n a n b ) (5 cd - n c n d ), n. 



2\ ' LIU ~ Lb - U / \~ L.U - -V I — * 



(7.17) 

Therefore at great distances we can write 



uv n=0 ri..r n 



l h$ ra (T,a) = -e — E A rsuv(n) E E n n- n r n — |yj + 0(l/r 2 ). 

(7.18) 
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Eq.(7.18) coincides with Eq.(3.34) of Ref.[21] 10 and the solution can be put in the form 



A hf^ rs (r, a) 



uv 



n) 



q\ uv (r- 






— * 

a 





+ 



+ (higher multipoles) + 0(1/ r 2 ). 



(7.19) 



Then Eq.(B13), i.e. q\ uv = \d 2 q uv \ TT ) leads to the standard quadrupole emission formula. 
Again one could evaluate the contribution of higher multipoles (mass octupole, momentum 
quadrupole,...) as in Ref.[21]. 

3. The Far Field of Time- Dependent Sources 

Let us look at the far field expression of our linearized 4- metric, given in Eqs.(6.12) in 
3-orthogonal gauges. 

By assuming matter sources with compact support and by using the multipolar expansion 
of Appendix B and Eq.(7.9), this equation can be rewritten in the following form in the far 
wave zone (r = |er|) 



A d T o 

!-7- 2 Z K '» 



— e 



M r + 



0, 



r 3 



A 



K, 



(i) 



A =—r{ M w c + J2^ 



AnG (ax j(i)) r ro , „, / , j \ 
_ 1_ \p T ando T (mass quadrupole) + 



M r = 

C r 
+ (higher multipoles)}, 



9(l)r 



= —eS r 



(mass quadrupole) + (higher multipoles)} 



(7.20) 



with the two asymptotic Poincare' charges M^c and j(i) given in Eqs.(4.21) and Eq.(4.23) 
respectively. Eqs.(B7) and (B12) have been used to get the result for 4 g(i) rr - The last 
term in A has been evaluated by omitting the electro-magnetic field and is negligible in 
the non-relativistic limit. The shift function M r gives the description of gravito- magnetism, 
Lense-Thirring effect included (see for instance chapter 6 of Ref.[2]), in the non-harmonic 
3-orthogonal gauges. 

We see that the results in the 3-orthogonal gauges are of the same type as in the standard 
harmonic gauges as can be seen by comparing Eq.(7.20) with Eqs. (13.30) and (13.32) of Ref. 
[27]. The only new terms are those involving the inertial gauge variable 3 /C(i) = -^ 3 K^. 



10 In Ref. [21] it is shown that for small velocities inside the source of gravitational waves the temporal 
derivatives of the stress tensor multipoles are negligible giving terms of order 0(v 2 /c 2 ). 
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Eqs.(7.20) are compatible with the (direction-independent) boundary conditions at spa- 
tial infinity for the 4-metric of our class of asymptotically flat space-times, given after Eqs. 
(2.7) of paper I (see also Eqs. (5.5) of Ref.[6]). 

C. PM Gravitational Waves with Asymptotic Background: the Energy Balance 

Having found the GW's of the HPM linearization of gravity, the next step is to check the 
energy balance: the energy emitted by matter in the form of GW's must be present in the 
gravitational field and there should be a back-reaction on matter. 

In the standard approach with compact sources there are various way to evaluate the 
energy balance: 

1) One can introduce the Landau-Lifschitz energy- momentum pseudo-tensor t^ L of the 
gravitational field and use the conservation law <9 M [— A g (T^ u + t^ L )]—0 to evaluate dE/dt 
in the far-field zone by using the quadrupolar approximation for GW's (see for instance 
Ref.[27]). This method is also used in the MPM + MPN approach of Refs.[24, 28], where 
the back-reaction (starting at the 2.5PN order) of the GW's on the source can be taken into 
account till 3.5PN. 

2) The same results are obtained with ADM Hamiltonian methods in Ref. [29] by taking 
a time-average of the work done by the quadrupole radiation-reaction force appearing in 
the Hamilton equations of the particles. Here it is emphasized the analogy with electro- 
dynamics due to the appearance of the analogue of the Schott term and it shown that there 
are no runaway solutions. 

3) Instead in Ref. [21] the coarse-graining method is used to find an effective energy- 
momentum tensor for the gravitational field at the second order, from which the increase in 
the energy of the gravitational field due to the emission of GW is evaluated. 

All these methods give the same result. The complications come from the problems of 
regularization of the gravitational self-force [16, 17] in the evaluation of the back reaction. 

Here we will show that we can recover the standard result without having the gravitational 
self-force, due to the Grassmann regularization {rjf mf = 0) of the gravitational self-energies, 
by using the conservation of the ADM energy (4.21). The effect will result from interference 
terms r\i rrii i]j rrij with i ^ j like it happens for the Larmor formula of the electro-magnetic 
case if Grassmann-valued electric charges are used to regularize the electro-magnetic self- 
energies [30]. For this calculation we consider only point particles as matter ignoring the 
electro-magnetic field, because we want to make a comparison with the treatments with 
compact sources (the calculations with the electro-magnetic field should add the assumption 
that such a field is localized in compact regions). 

By making some integrations by parts (assumed valid with our boundary conditions), the 
ADM energy (4.21) has the form (M% V) is given in Eq.(3.12)) 
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- E ADM = j d*a (M$P + M%P + ^ M<™ 1 M<™ 



c 



+ 



+ 



c 3 



16nG 

ab 



[dr R- a M- al T R- b + Y, d a R- a M- ai d a R,] ) (r, a) + 0(( 3 ). 

a 

(7.21) 



Since we have: 

1) a* = E. 7a.it> (E r if> = o) ; 

2) M ab = EaS To-Tw^ from Eq.(6.4); 

3) E b M«*r« = -§ Xi)aa from Eq.(6.8); 

4) 5 rs r^ = -§ ( 4 /i(i) rs - |5 rs E„ 4 fyi)™) fro m 4 ^(i)rs = -2e<5 rs (r^ 1) + 2<f>^) with 

2< ^(i) = _ l E„ 4/l (i)^; 

5) 5 rs r« = -f ( 4 /^ rs + | (9 r e (1)s + d s e (1)r ) + («9 r <9 S - § <5 rs A) A (1) ) from Eq.(6.5); 

6) Er 4 ^ = E^ r 4 ^ s = 0; 

the term bilinear in the gradients of R a in Eq.(7.21) can be manipulated in such a way that, 
after an integration by parts, the final form of Eq.(7.21) becomes 

i Eadm = / *a (MfP + M%P + ^ <"> ^ <"> " 

- E <>? [| 3 %> - ^ { «>:> - \ E a -§ *S?) 

c d 

+ iSG E + E ( a '" h mr,) 2 ] )(r,3) + O(C') = 



+ 



rs 

^ 1 | d 3 ap E{1+2) (r,a) + 0(C 3 ), (7.22) 

where in the last line we introduced the (coordinate dependent) density Pe(i+2)(t,<t) of the 
ADM energy E ADM till the order 0(( 2 ). This energy density is the sum of three terms: a 
matter term P^i+2)^( r ' a radiation term /0^n^. 2 )( r ' ^) (^ ne GW) and an interaction term 
Pb(i+2) ( r ' ^0 ( tne interaction of GW's with matter: it controls both the emission and the 
back- reaction). Therefore we have 

/ -»\ (matter) i ->\ . (rad) / ->\ , (int) / ->\ 

PE(l+2) (T, 0") = #5(1+2) ( T ' CT ) + Pi5(l+2) ( T - a ) + Pj5(l+2) ( T ' CT ) > 
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P { e^(t, a) = M {1) (t, a) c = ^ 5 3 (<t - ^(r)) Vi c sj mf c 2 + /? 2 (r), 



/27TG 



,(c/v) 



,([/V) 

(i) 



" ; v 1 / M (uv) 1 dcdd M (uv) 



M) 



-E 



G 



E ^ / rf 



r £" + 5E| r«)(r, 4 (r))- 

d 



ffi(T) -a\ \B-rjj(j) 



6 jV« 



+ 



G 



8tt c 2 



E ^( r ) E Kid ^ E ^ / d3(j t^tt — =r ^ ^ r= — ^7~TT ~~ 



G ^ V m l c 2 + «? 0") a/ ™ 2 c 2 + k 2 (t) 

— E — 



2 c 



(7.23) 



Let us divide the 3-space S T in two regions by means of a sphere 5 of big radius R» l c : 
a) an inner region V^ nner ) with a compact sub-region of linear dimension l c containing all 
the particles (and the electro- magnetic field if we would add it); b) an asymptotic far region 
V(/ ar ). Let rf = cr r /\a\ be a unit 3- vector. 



Since we have Eadm = E}^m + E 



ADM ADM ' an d P£(l+2)H r > °"J keV /ar — Ps(l+2) a )\v€V f 



(matter) 



(int) 



0, we get ^ = / V(/ap) rf 3 ap^| 2) (r,a). 

Since Eadm is a constant, we have dE ^ M = so that we get 
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u ^ 'ADM u ^ADM 



dr dr 



c 4 



327TG 



C 



32ttG 

.4 



C 



32ttG 



-/ d 3 ad T p { ^l 2) (T,a) = 

JV (far) 

d*a [Or % T )rs 81 % T )rs + £ d c % T )rs d T d c ^ rs ] (r, a) = 

r) rs c 

[ [E^ d c % T )rs ) + T ^ rs □ % T )rs ] (T, <?) = 

/ ^EE ^ (^ 4 ^W«5 c 4 ^ rs )(r,a) + 0(l/r 2 ), (7.24) 



'(/<"■) rs c 

where we have done an integration by parts and used Eqs.(6.11) (their second member is 

(i) 



zero because the non- local TT quantity T^ T)rs {T, a) is of order 0(l/r 2 ) (r = \a\) in V (/ar) ). 



If in the far region we use the solution (7.19), we have d c 4 hj^ rs (r,(x) = 
-n c d T %l s (r,a) + 0(iy) and £ c n c d c % T )rs (r, a) = -d T % T )rs ^, <?) + 0(l/r 2 ) be- 
cause r=|cr|>i?>>/ c . Therefore we get (the sphere S is parametrized with the angles 9 
and ip) 

J Thinner) j TTi^t/ar) 



_ u ADM 



c 4 



327r G ^ , 

= 3^G R2 I d ( cosd ) d ^T,{ d r%ls)\r,a) + 0(l/r 2 )^ 
J s rs 

T i<i: 6 / d(cos9)d<p V A^n^g^r-i^g^r-i?) + 0(l/r 2 ) = 



2 



2/T C 

= ITS E ^ ^ S|rT ( r _ ^ ^ ^' Tr ( r _ ^ + °(V^)- (7-25) 

rs 

where in the last line we used Eq.(B13) and J s d(cos 6) dip A rsuv (n) = || (11 S ru 5 SV — 
A5 rs S uv + 5 rv 5 SU ) (see Eq.(3.74) of Ref.[21]). The change of sign is due to the fact that 
the unit vector n c is minus the normal to the sphere S. 

Therefore Eq.(7.25) reproduces the standard result for the total radiated power also 
named the total gravitational luminosity C of the source (see for instance Eqs. (1.153) and 
(3.75) of Ref. [21]) 11 . 



j g (inner) / rs \ TT \Z 

Since r = ct, Eq.(7.25) can be rewritten as — — = -^r ( q at - 6 — J . However, our mass quadrupole 
qrs\TT - 1S e q ua i ^ cQ rs , where Q rs is the mass quadrupole of Ref. [21]. As a consequence, we have 
(^P) 2 as in Ref. [21]. 



ADM Ct 



dt 5 c 5 
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See Appendix C for the evaluation of the balance equations for the 3-momentum and the 
angular momentum by using the conservation of the corresponding ADM generators (4.22) 
and (4.23). 

D. PM Gravitational Waves with Asymptotic Background: Detection with the 
Geodetic Deviation Equation 

In Ref.[21] there is a full discussion of the detectors of gravitational waves and of the 
reference frames for the observers looking for them. After discussing the observers using 
local inertial frames (with Riemann normal coordinates) or freely falling frames (with Fermi 
normal coordinates), there is a discussion of observers using a TT frame (where particles 
at rest before the arrival of the gravitational wave remain at rest after their arrival) and 
of their use of the geodesic equation and of the geodetic deviation equation. Then there is 
a discussion of the proper detector frame used by experimentalists (in it locally one uses 
an Euclidean Newtonian 3-space): now the solution of the geodetic deviation equation, 
giving the coordinate displacement £ (t) of a test mass of mass m induced by a gravitational 
wave (zero in the TT frame), can be put in the form m d = F l with the Newton force 

d 2 4 h TT ■ 

F % = y dt 2 i 3 ■ However, it is possible to give a coordinate-independent description of the 
effect of a gravitational wave on a test mass by using a notion of proper distance between 
two nearby geodesies, see Refs.[20, 22]. 

All these approaches could be reproduced with our formalism. In this Subsection we give 
first a discussion of geodesies in our generalized TT gauge as our alternative to the TT frame 
and then we discuss the proper distance and the geodesic deviation equation. 



1. The Geodetic Equation for Test Particles 

Let us consider the behavior of a test particle in presence of a gravitational wave in the 
distant wave zone far away from the sources. To this end let ignore the sources by putting 
M.^^ = 0, M-^^J = 0, Tffj = (they are assumed to have compact support and to give 
vanishing terms E a = in Eqs.(6.4) evaluated in the far zone) and consider a gravitational 
wave solution of the wave equation □ hj^ aa (r, a) — 0, see Eq.(6.8). 

In this case the 4-metric is given in Eq.(6.12) with all the matter terms omitted (so that 
we have (f)^ = ^ c ^ r^) and with 3 /C(i) = (here we assume that this 3-orthogonal 
gauge is relevant near the detector). 

By using the embedding z M (r, a) = e p A a A discussed in the Introduction, we describe the 
world-line x^(s) = e^cr A (s) = e^r(s) + e^ff{s) of the test particle as a time-like geodesic 

with parameter s ( d + 4 T^ c (a E (s)) da d ^ dc = 0). If we choose for s the proper 

time of the test particle, we have s = s(t) = \J e A g(i) TT {r(s) 1 r)(s)) and £ M (s(r)) = x^{t) = 
e\ V A (r) = r + & r rf (r ) . We have 

- = l + n (1) (r,^(r)) + 0(C 2 ). (7.26) 
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In the weak field approximation the geodesic equation becomes the following equation 
for fj(r) 



d 2 n r (r) 
dr 2 



= -( 4 r[ 1)rr + 2 4 q 1)ru ^+ 4 q 1)tt ^v) + 



+ ( 4 r[ 1)rr + 2 4 r[ 1)ru ^ + 4 r[ 1)w ,i)V) v r + \o{C). 

Let us consider the following Christoffel symbol 



(7.27) 



6 i (l)TT 



<9n ( i) , 9n ( i )(r) + I ^2) 



+ 



<9r 



— arw - — V Srw 

A a 9 A A c 



+ ^o(c 2 ) 



= a^ 1) + o(0, 



(7.28) 



where we used the function defined in Eq.(6.13). 

As a consequence the geodetic equation can be written in the form 



dr 2 



( l > \ -u I 4-pr ^, I y.wy.^ I 



(i) 



+ ( 4 r[ 1)rr + 2% )TU v u + % )uv r, u V v ) V r + T 0(C 2 )- (7.29) 

A special set of geodesies solution of Eq.(7.29) is implicitly defined by the condition 

r/ r (r) - ^ (1) (r,7f (r)) = constant. (7.30) 



In this case we have 



dr da s 

so that in the weak field approximation we get 



(7.31) 



(7.32) 



As a consequence the right side of Eq.(7.29) is zero modulo terms of order j-0(( 2 ): 



dr 2 



^M-^(M r W)]=o + |o(( 2 ). 



Moreover in the weak field approximation the 4-velocity of the geodesic (7.30) turns out 
to be 



a da A (s) dn A {r) dr(s) ( r dVfa 

U = ; = ; ; = ( U = 1 — 71(1); U ~ 



ds 



dr ds 



dr 



(7.33) 
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Let us remark that, since we have rf = + 0(( 2 ) = 0(() along these geodesies, then 
we get jfc = §J + §^fj s = §^ + \ 0(C 2 ) for every quantity A(t, a u ) of order 0(C). 

The solution (7.30) selects a special family of geodesies whose meaning can be clarified 
by remembering the coordinate transformation (6.13) leading to a generalized TT gauge, 
whose associated 4-metric is given in Eq.(6.15) with the matter terms omitted. 

The new Christoffel symbols 4 f A BC = g£ ( 4 T D EF §fi + ^gcr) imply 

r ( r 1)rr = + io(C 2 ). (7.34) 

This equation is the sufficient and necessary condition to get the result that the coordinate 
lines cr r = constant be geodetic. This consequence of the solution (7.26) is in accord with 
the choice of constant spatial coordinates for the usual TT gauge done in Ref. [20], pp. 
13-16. 



2. Detection of Gravitational Waves 

As shown in Ref. [22] the main observable for the detection of gravitational waves is the 
proper distance between two nearby geodesies. 

The geodetic deviation S a (t) is the infinitesimal 4- vector orthogonal to the 4- velocity 
u a {t) of the reference geodesic: ua(t) £ a {j) = 0. If n A {r) = (r, n r {r)) is the reference 
geodesic and n A (T)+£ A (r) the nearby one, the proper distance between them is the invariant 

D(t) = ^g AB (T,n-(T))£ A (r)£B(r). (7.35) 

Following Ref. [20] and consistently with the weak field approximation we assume 
D(t) « L, namely that the proper distance is less of the wavelength of the gravitational 
wave to be detected. 

Given the geodetic deviation equation for £ a (t) 

(u a (t) V a ) {u b {t) Vb) 8 C {r) = - 4 R C 'efg{t] D {t)) u E {r) £ F {r) u G (r), (7.36) 
we get the following equation for the proper distance (a scalar quantity) 



(u A (r)V A fD(r) = ^ D(s(r)) 



- 4 Rabcd(v E (t))u b (t)u d (t) 



£ A {r)£ c {r) 



D(r) ■ 

(7.37) 



In the weak field approximation Eq.(7.37) becomes 
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(l - 2n (1) (^ M ) 4^ - D (r) = 

= - 4 B (1) „.„(r,^"(r)) +\0(C\ (7.38) 

with the following expression for the Riemann tensor implied by Eq.(6.12) 

4o _ 92 n W 19 ( dn {l){r) dn (1)(r) \ d 2 , (1) , 1 2 _ 



1 /) 2 fc TT _L_ j_ 

2 °V ft (l)rs + ^2 



d 2 r hJl s + -0{C). (7.39) 



By using the result 7^ = §7 + x 0(C 2 ), valid for every quantity A(t, a u ) of order 0((), 
we get the following final form of the equation for the proper distance 

(l _ 2 ^ r)) )*^ ) _-S Sa M- D(T) . 

+5 £^.™> £ ^ + !<>«'>• (7.40, 
To study Eq.(7.40) we follow the method of Ref.[20]. 

The structure of Eqs. (7.36), (7.37) or (7.40) suggests that the geodetic deviation £ A (r) 
may be parametrized as the sum of a constant deviation £ A plus a small corrective term 
determined by the weak gravitational field. To this end let us introduce an arbitrary constant 
direction with the unit constant 4-vector n A , 4 gAB nAnB = 1, so that the constant part of 
the deviation is given by £ A = D Q n A with D Q = const.. If 5£ a (t), with 5£ D ^ = 0(C), is 
the small corrective term, then the geodetic deviation is parametrized in the following form 

£ A (r)= D (n A + 6 -^l + 0(C 2 ))- (7.41) 



o 



The induced parametrization of the proper distance is 



D(t) = D (l+ ^ + 0(C 2 ) ) , (7.42) 



D 



with ^ = Vmb^^ = 0(C). 

With these parametrizations Eq. (7.40) takes the form 
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allowing to find the solution for — ^ at the lowest order. 

As a consequence we get the following expression for the proper distance 

D(r) = D Q + + 1 - h^ rs (r, V u (r)) rf n s + 0(( 2 ). (7.44) 

The choice of the direction n r allows to discuss the effects of the polarization of the 
gravitational wave on the detectors as it is done in Ref.[20, 22]. 
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VIII. CONCLUSIONS 



We have denned a consistent HPM linearization of ADM tetrad gravity in the York canon- 
ical basis in the family of non-harmonic 3-orthogonal Schwinger time gauges parametrized 
by the numerical value 3 Km(r, a) of the York time, the inertial gauge variable describing 
the general relativistic remnant of the freedom in clock synchronization. The non- Euclidean 
instantaneous 3-spaces S T (a first order deformation of the Euclidean 3-spaces of inertial 
Minkowski frames) are dynamically determined except for for the value of the trace 3 K^ of 
their extrinsic curvature tensor. The 4-metric has an asymptotic Minkowski background at 
spatial infinity. PN expansions are avoided by introducing a ultraviolet cutoff. 

The PM solutions, 0(i)(t, a) — 1 + 6 0(i)(r, <t) (the 3- volume element), 1 + n^(r, a) (the 
lapse function), f2,(i)( r )(r, a) (the shift functions), o r (i)( r )( s ) \v^s{t, <?) (the non-diagonal ele- 
ments of the shear of the Eulerian observers), of the constraints and of the Hamilton equa- 
tions implying the preservation in r of the gauge-fixings, have both an instantaneous action- 
at-a-distance part depending on the instantaneous value of the matter energy-momentum 
tensor and a part depending on the tidal variables. At the PM level all previous quan- 
tities do not depend on the York time 3 Km (r, a) but on the spatially non-local function 
3 /C(i)(t, cr) = ^ 3 fT(i)(r, a) (it can be named non-local York time). 

In these non-harmonic gauges two functions M a i{(?) Ra{T, (t), with the operator 
M a i(a) containing only spatial derivatives, of the tidal variables R a (r, a) satisfy a wave 
equation with the flat asymptotic d'Alambertian. It has been shown that the operator 
M a i(<j) contains the information for determining the TT part of the 3-metric on E T . By 
using a no-incoming radiation condition with respect to the flat asymptotic metric, we get 
a retarded solution for the tidal variables in terms of the matter energy-momentum tensor 
which describes PM TT GW's with asymptotic background propagating at the velocity of 
light c 12 . A multipolar expansion of the energy- momentum tensor in terms of Dixon mul- 
tipoles allows to get the standard quadrupole emission formula as the dominant part. Also 
the rate of variation of the energy is correct due to the Grassmann regularization 13 of the 
self-energies of the point particles (no gravitational self-energy) due to the existence of the 
conserved ADM energy. Thereforall the main properties of GW's are reproduced in our 
HPM approach. Only the hereditary and memory tails (coming from the matching of MPM 
and MPN solutions as said in Appendix A) are missing at this order: to study them we have 



See Ref.[31] for the problem of the velocity of propagation of the gravitational field in general relativity 
and in bimetric theories. 

In the electro-magnetic case [15] the Grassmann- valued electric charges allow to find the effective potential 
(Coulomb plus Darwin) corresponding to the one-photon exchange Feynman diagram. Therefore the 
problem of electro-magnetic self-energies is pushed to the level of loop diagrams and becomes part of the 
problem of renormalization of QED. At the classical level the Grassmann regularization allows to make 
sense of the classical equations of motion (not well defined due to essential singularities at the charge 
location): the replacement of the Grassmann- valued electric charges with the standard electric charge is 
equivalent to a classical renormalization of scalar charged particles. With gravity the same mechanism 
is at work, except that we do not yet have an accepted renormalizible theory of quantum gravity. Our 
procedure is a classical regularization and in some sense we are identifying the effective potential connected 
to the one graviton exchange. 



80 



to go to higher orders in the HPM expansion (see Section IIIB), which should correspond 
to the MPM solutin in the far wave zone. 

All these results lead to a PM solution (modulo the choice of the inertial gauge variable 
3 /C(i) = ^ 3 fT(i)) for the gravitational field and identify a class of PM Einstein space-times. 

These results are gauge dependent because the York time 3 K^ is an inertial gauge vari- 
able. Therefore we have to face the gauge problem in general relativity. The gauge freedom 
of space-time 4-diffeomorphisms implies that a gauge choice is equivalent to the choice of a 
set of 4-coordinates in the atlas of the space-time 4-manifold. 

The standard approach to the gauge problem is to try to describe physical properties 
in terms of gauge invariant quantities, i.e. in terms of 4-scalars. At the Hamiltonian level, 
where the Hamiltonian gauge group (whose generators are the first-class constraints) is equiv- 
alent to the 4-diffeomorphism group of space-time only on-shell (i.e. on the solutions of Ein- 
stein equations; see for instance Refs.[5]), the Hamiltonian gauge- invariant (off-shell) quan- 
tities are the Dirac observables (DO), which have zero Poisson bracket with the constraints. 
If we would know the solution of the super-Hamiltonian and super-momentum constraints 
in the form = cf) — F[9\ ni, Rg,, n s , matter] 0, nf^ = itf^ — Gi[9\ tt i,R a ,Ila, matter] ps 0, 
we could look for a canonical transformation from the York canonical basis to a Shanmu- 
gadhasan basis containing n, n( r ), 9\ Tt\ • ~ 0, ~ 0, n^, R a , H a , DO(, matter ). In this basis 9 l 

and ir^, n, n^, would be the (primary and secondary) inertial gauge variables and R a , Il a , 
DO (matter) the DO's (non-local quantities as functions of the original variables). However 
no-one is able to find such a basis. Also the most recent works of Refs. [32] contain existence 
proofs but no workable algorithm for finding the Dirac observables of the gravitational field. 
Moreover it is not clear how many of these DO are 4-scalars. Hopefully the 3-scalar tidal 
DO's R a) Tl a , may be replaced with two pairs of 4-scalar DO's connected with the eigenval- 
ues of the Weyl tensor [33] by means of a canonical transformation (this conjecture is under 
investigation). For the transverse electro-magnetic field one expects that the final DO's can 
be chosen as tetrad-dependent 4-scalars by using the Newman- Penrose formalism [33]. 

On the other side at the experimental level the description of baryon matter is intrin- 
sically coordinate-dependent, namely is connected with the conventions used by physicists, 
engineers and astronomers for the modeling of space-time. 

A) The description of satellites around the Earth is done by means of NASA coordinates 
[34] either in ITRS (frame fixed on the Earth surface) or in GCRS (frame centered on the 
Earth center) (see Ref.[35]). 

B) The description of planets and other objects in the Solar System uses BCRS (a quasi- 
inertial Minkowski frame, if perturbations from the Milky Way are ignored 14 ), centered in 
the barycenter of the Solar System] and ephemerides (see Ref. [35]). 

C) In astronomy the positions of stars and galaxies are determined from the data (lumi- 
nosity, light spectrum, angles) on the sky as living in a 4-dimensional nearly-Galilei space- 
time with the celestial ICRS [36] frame considered as a " quasi- inertial frame" (all galactic 



Essentially it is denned as a quasi-inertial system, non-rotating with respect to some selected fixed stars, 
in Minkowski space-time with nearly-Euclidean Newton 3-spaces. The qualification quasi-inertial is in- 
troduced to take into account general relativity, where inertial frames exist only locally. It can also be 
considered as a PM space-time with 3-spaces having a very small extrinsic curvature. 
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dynamics is Newtonian gravity), in accord with the assumed validity of the cosmological 
and Copernican principles. Namely one assumes a homogeneous and isotropic cosmological 
Friedmann- Robertson - Walker solution of Einstein equations (the standard ACDM cosmo- 
logical model). In it the constant intrinsic 3-curvature of instantaneous 3-spaces is nearly 
zero as implied by the CMB data[37], so that Euclidean 3-spaces (and Newtonian gravity) 
can be used. However, to reconcile all the data with this 4-dimensional reconstruction one 
must postulate the existence of dark matter and dark energy as the dominant components 
of the classical universe after the recombination 3-surface! 

As a consequence of the dependence on coordinates of the description of matter, our 
proposal for solving the gauge problem in our Hamiltonian framework with non-Euclidean 
3-spaces is to choose a gauge (i.e. a 4-coordinate system) in non-modified Einstein gravity 
which is in agreement with the observational conventions in astronomy. Since ICRS has 
diagonal 3-metric, our 3-orthogonal gauges are a good choice. We are left with the inertial 
gauge variable 3 /C(i) = ^ 3 -^(i) (not existing in Newtonian gravity). The suggestion is to try 
to fix 3 /C(i) in such a way to eliminate dark matter as much as possible, by reinterpreting it 
as a relativistic inertial effect induced by the shift from Euclidean 3-spaces to non-Euclidean 
ones (independently from cosmological assumptions). As a consequence, ICRS should be 
reformulated not as a quasi-ineriial reference frame in Galilei space-time, but as a reference 
frame in a PM space-time with 3 K^ deduced from the data connected to dark matter: as 
a consequence what is called dark matter would be an indicator of the non-Euclidean nature 
of 3-spaces as 3-submanifolds of space-time (extrinsic curvature effect), whose internal 3- 
curvature can be very small if it is induced by GW's. Then automatically BCRS would be its 
quasi-Minkowskian approximation (quasi- inertial reference frame in Minkowski space-time) 
for the Solar System. This point of view could also be useful for the ESA GAIA mission 
(cartography of the Milky Way) [38] and for the possible anomalies inside the Solar System 
[39]. 

As a consequence in the third paper [18] we will study the PN expansion of our HPM 
solution for a system of point particles without electro-magnetic field. There we will de- 
termine the explicit dependence of the equations of motion of the particles, of the proper 
time of time-like observers, of the time-like and null geodesies, of the redshift of light and 
of the luminosity distance upon the time- and spatial-gradients of the non-local York time 
3 /C(i) = 4^ 3 if(i)) in the PM space-time. Then we will study the slow motion limit making a 
PN expansion at all ^PN orders with a detailed study of the order 0.5PN. 

We will show that the main observational supports for the existence of dark matter 
(rotation curves of galaxies, mass of galaxies from the virial theorem and from gravitational 
lensing; see for instance Refs.[40]) can be translated in restrictions on the non-local York 
time 3 /C(i)(r, a) at the location of galaxies. If we could find a global phenomenological 

parametrization 3 /C^ en ' ) (r, a) over all the 3-universe (in our PM space-time and with some 
average for the r-dependence), we could get a phenomenological determination of the York 
time 3 IC { ^ en \r,a) = A 3 JC^ en \r,a) to be used as an observational clock synchronization 
to be used to define the 3-spaces of a PM ICRS. 
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Appendix A: The Standard Approach to Gravitational Waves by using Einstein's 
Equations in the Family of 4-Harmonic Gauges 

The standard description of GW's (see for instance Ref.[21]) is done in the family of 
4-harmonic gauges after a decomposition of the 4-metric 4 g iiV of asymptotically flat space- 
times in a flat Minkowski background 4 i] fJiU plus a small perturbation, 4 g^ v = 4 r\ llv + 4 h fiu , 
l^h^l << 1. The linearized Einstein equations in harmonic gauges have the form \D 4 h^ u = 
-e^T^, d v4 W v = ( 4 V = 4 V - \ \»% 4 h = VXo e = ± according to the 
signature convention for the 4-metric), where T^ v is the matter energy- momentum tensor 15 
satisfying d u = at this order. 

It turns out that with an allowed coordinate transformation inside the family of harmonic 
gauges, i.e. 4 h^ v H> 4 /i^ = 4 h^ v - (<9^ + d u ^ - 4 ^<9 P £ P ), one can identify a transverse- 
traceless (TT) harmonic gauge in which the only remnant of the gravitational field are the 
two polarizations of GW's. 

This weak field description is assumed valid in the wave zone far away from the matter 
sources, which are assumed in slow motion (Post-Newtonian (PN) approximation v « c) 
and with a small self-gravity 16 (if d is the linear dimension of the source with mass M 
and R M = its Schwarzschild radius, one has ^f- « (^) 2 = e << 1, i.e. d » R M ] 

e n/2 = (2)n is t he | p N orc l e r). 

In this way one is led to the interpretation of a (spin 2) GW propagating in the Euclidean 
3-spaces of an inertial frame in Minkowski space-time. 

Since the equivalence principle forbids the existence of global inertial frames in Einstein 
space-times, the next step is to replace the above decomposition with one with respect to a 
curved background 4 g fiu , i.e. 4 g lxv = 4 g fiu + 4 h ( jfu , such that the GW's are only ripples on this 
background 17 . With this formalism and a coarse-grained description 18 one can evaluate 
the energy loss associated to the emission of GW's and take into account the back-reaction 
on the background which is therefore modified. 

However the most advanced description of GW's is done with the Damour-Blanchet ap- 
proach [24] or with the equivalent DIRE approach of Ref.[28] (see chapter 5 of Ref.[21] 
for an overall review; see also Ref.[20]). In these approaches one uses the decomposition 
4 fiV = + A h^ v and writes a so-called relaxed form of Einstein equations for the quantity 
4 h^ = \J\ 4 g\ 4 g^ v — 4 rj^ v , namely D 4 h Miy = e t^ u , where the effective energy-momentum 



OurT M „ is ±T M „ of Ref.[21]. 

See Ref.[21] for a review of the problem of self-gravity of extended compact objects, in particular binary 
systems, and of the effacement principle of the internal structure which becomes relevant only at the order 
5PN. 

As shown in chapter 1 of Rcf.[21] this decomposition makes sense if: A) cither i g flL , has a scale of spatial 
variation Lb and the wavelength of the GW is A << Lb (and \ 4 h^\ « X/Lb « 1); B) or 4 g flv contains 
only frequencies less than vb and the frequency of the GW is v » vb (this case is the more relevant for 
detectors). Moreover, a detector of dimension Lp will react only to GW's with A >> Lb- For GW's with 
frequency 10 -4 — 1 Hz (to be detected by LISA) the wave-length is of order 10 6 — 10 10 cm (A v = c). If 
the frequency is 1 — 10 4 Hz (to be detected by LIGO, VIRGO), the wave-length is of order 10 10 — 10 14 cm. 
Either a spatial average on many wavelengths A of the GW or a temporal average on several periods 1/v 
of the GW (the method used in the detectors). 
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tensor is r"" = | 4 #| + ^A^Ch) with A^( 4 h) = ^ | 4 #|C + 9 a 4 h^^ 4 h^ a - 
4 h a/3 <9 a 9g 4 h MI/ . t^ L is the Landau-Lifschitz energy- momentum pseudo-tensor of the gravi- 
tational field, which satisfies the ordinary conservation law d v \\/\ 4 g\ (T^ + tj^ L )] = due 
to Einstein equations. The ordinary equations of motion for the matter (i.e. 4 V„ T^ v = 0) 
are obtained by restricting the solutions of the relaxed Einstein equations to the harmonic 
gauges by requiring d v 4 h^ u = 0. 

Since the matter is supposed to have compact support of size d and to be in slow motion 
(y^e = - c Rj \J~^f- « 1 and typically with a wavelength of GW's satisfying A >> d), the 
solution of the relaxed Einstein equations is obtained in three steps: 

A) In the far wave zone (r >> A and d < r < oo), where T^ v = 0, one makes a Post- 
Minkowskian (PM) expansion 4 h^ = Yl^=i G n 4 h^ and uses the restriction d u 4 h^ = 
to harmonic gauges. The iterative PM solution (including a homogeneous solution of the 
wave equation) is expressed in terms of retarded (symmetric trace free or STF) radiative 
multipoles of the gravitational field. Only a finite number of multipoles in 4 h± u are taken 
into account to avoid the problem of the Green function of the wave operator, which would 
require the knowledge of 4 h^ also in the near region where the PM expansion does not hold. 
With this regularization a multipolar PM (MPM) solution is found with the property that 
by making a PN expansion one finds 4 h°° = 0(l/c 2n ), 4 h™ = 0(l/c 2n+1 ), 4 hjf = 0(l/c 2n ). 

B) In the near zone (r << A; the exterior near zone is d < r << A) one makes a PN 
expansion of 4 h fiu with 4 h 00 = J2n=i 4h oo\ 4/i oi = J2n=i 4h oi l+1 \ 4}l ij = E^=i with 

= 0((f ) n ) (v/c « a/ Rm I d) and of the energy-momentum tensor T 00 = J2n=o ( 2 ™)T 00 , 
rpOt _ ^2™ =l (2n+i)yo«^ rptj _ Y^ =1 ( 2n )T u . Since the Newtonian approximation corresponds 
to 4 /4o\ 4 hoi — 0, 4 hij = 0, the 1PN order contains 4 h$ , 4 h^\ 4 h^\ the 2PN order contains 

4 h^Q , 4 hfi , 4 h^\ and so on. Beyond some order one finds divergences connected to the 
inversion of the Laplacian operator which require the introduction of a regularization of 
the Poisson integrals. This is due to the fact that it is not possible to rebuild a retarded 
solution from its expansion for small retardation without going outside the near zone (i.e. 
there are divergencies for r — > oo). To avoid these problems one introduces a multipolar PN 
(MPN) expansion and uses it in the relaxed Einstein equations with 4 \^ v = X]^=2 ^ 4 ^(n)> 
T tiu _ Y^ = _ 2 j™ T (n)- By using a regularization prescription one finds retarded solutions 
(regular at r = 0) in terms of STF matter multipoles. Also a homogeneous solution of the 
wave equation is needed: its regularity at r = requires that it is a half retarded minus half 
advanced solution. Finally to get the equations of motion of matter and their PN expansion 
till order 3.5PN one has to impose the harmonic gauge condition and to take into account 



(2n+l) 
00 



the back-reaction from the emission of GW's: this introduces the lacking terms 4 h, 
4 4f \ 4 hf j n+1 \ and becomes relevant at the 2.5PN order [0((f) 5 )]. 

C) Then one matches the MPM and the MPN solutions in the overlap of the near and 
far zones: this allows to express the radiative multipoles in terms of the matter ones. Now 
one can study the limit at future null infinity (r — > oo with u — r — t/c fixed) to test the 
nature of GW's. At higher orders hereditary terms (tails starting from 1.5PN [0((^) 3 )] and 
non-linear (Christodoulou) memory starting from 3PN (see Ref. [41] for a review) appear, 
showing that GW's propagate not only on the flat light-cone but also inside it (i.e. with all 
possible speeds < v < c): there is an instantaneous wavefront followed by a tail traveling at 
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lower speed (it arrives later and then fades away) and a persistent zero-frequency non-linear 
memory. 

Today there is control on the solution and on the matter equations of motion till order 
3.5PN [0((^) 7 )] (for binaries see the review in chapter 4 of Ref.[21]) and well established 
connections with numerical relativity (see the review in Ref.[42]) especially for the binary 
black hole problem (see the review in Ref.[43]). 

In Refs.[29, 44] there is a Hamiltonian approach starting from a PN expansion of the 
ADM formalism in suitable non-harmonic gauges (generalized isotropic ones with 3 K = 0), 
which allow to recover the previous results in harmonic gauges till the order 3PN. 

However in this formulation GW's propagate in the background Euclidean 3-space implied 
by the decomposition A g fiu = A i] fiu + 4 h fll/ and it is not clear how to visualize them in the non- 
Euclidean instantaneous 3-spaces of the global non-inertial frames implied by the equivalence 
principle. 
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Appendix B: The Multipole Moments of the Matter Energy-Momentum Tensor 



In Ref.[23] there is a study of the relativistic Dixon multipoles [45] 19 of the energy- 
momentum tensor of relativistic matter systems (for instance point particles plus the electro- 
magnetic field) in the rest-frame instant form of dynamics in Minkowski space-time. Since 
in HPM gravity in our asymptotically flat space-times and in its non-inertial rest-frame 
developed in this paper we have only small deviations from such a scheme, we can apply 
this formalism (without using the general relativistic Dixon multipoles [48]) to the energy- 
momentum tensor given in Eqs. (3.12), which gives rise to the ten internal Poincare' gener- 
ators given in Eqs. (4.21) - (4.37) at lowest order 0(C)- We refer to chapter 3 of Ref. [21] for 
a review of the standard types of multipoles used for the study of gravitational waves from 
compact sources. 



Moreover at the lowest order we have d A Tff = 0, d A (Tff a c - Tff a B ) = 0, with Tfi = 
fj\ Tfi = MfJ r \ see after Eq.(3.12). 



Let w^ij) = 2 M (r, ff(r)) be the time-like world-line of an arbitrary center of motion of 
the particle system. On the instantaneous 3-space E r we can define the following Dixon 
multipoles of the matter energy- momentum tensor T£ b (t, a) with respect to the point t)(t) 

q r,..r n \ ab (t) = J Sa{a^-^{r)).. .(<^ - if»(r)) 7^(r,*), (Bl) 

and the following multipolar expansion of T^ B 



n=0 ri...r„ 

T^ T)ab (r, a) = ^)) (t, a), (B2) 



where Eq.(6.9) has to be used for the TT tensor. 

To connect them to the standard Dixon multipoles of the space-time energy-momentum 
tensor T txv (x = z(r,a)) = z^z^T^ B ~^ (r, <r) + 0(( 2 ) we must use the adapted embedding 
z^(r,a) discussed in the Introduction. 

By using Eqs.(4.21)-(4.37), the relevant multipoles are: 

la) the mass monopole 

q lTT = J d 3 aM^ V) (r,a) = M (1) c; (B3) 
lb) the mass dipole (Eq.(4.37) is used) 



19 As shown in this paper, strictly speaking the multipolar expansion makes sense only if the energy- 
momentum tensor is an analytic function of the 3-coordinates. However see Ref. [46] for the relaxation of 
this condition. 
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9 '|^/ A (,-, W ,<>M,. j;; -m (11c ,' W!! 

« -M (1) c^(r); (B4) 
lc) the mass quadrupole (Eq.(4.37) is used) 

q rs]TT = J d 3 a(a r -rf{r)){a s - V s (r)) M^ V \r,a) » 

« J d 3 aa r a s M^ ) V \r : a) + M {1) crj r (T)rj s (T); (B5) 

2a) the momentum monopole (Eq.(4.22) is used) 

g'- = | d 3 aM^ V r \r,a)=p r (l) ^0; (B6) 
2b) the momentum dipole 

d rWs = J d 3 a(a r - V r (r))M^ V s \r,a) « | ^^'Jfr,?), (B7) 

whose antisymmetric part is g r ' rs — g s ' rr = —2j^ due to Eq.(4.23) [as a consequence we 
have q r I rs = \ (q r I TS + q s I Tr ) - ; 

3) the stress tensor monopole 

q lrs = J d 3 aT$(T,a) = q\ sr . (B8) 

If we choose as center of motion of the mass distribution the center of energy w^(t), we 
must put equal to zero the mass dipole and this implies ^ r (r) «s (namely w^(t) = z^{t, 0) 
coincides with the origin of 3-coordinates) 

rs \ TT = J d 3 aa r a s M^p(r, a). (B9) 



q 



In this case the non-zero lowest multipoles are the mass monopole Mm c, the mass 
quadrupole g rs l TT ; the momentum dipole g r l TS and the stress tensor monopole g' rs . 

The lowest order conservation law 8aT^ = gives d T M^^ = — J2 r d r M.^P and 

drM^y^ = — Yls ^s^(T)' see Eqs.(3.13). By integrating over the whole 3-space with the 
matter density having compact support (or suitable fall-off at spatial infinity) we get 
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d T M (1) = + 0(C 2 ), d T pft = + 0(C 2 ). (BIO) 
The conservation law d A (T^ cr c - T^P a B ) = implies 

drift =0 + 0(C 2 ), drjft = + 0(C 2 ). (Bll) 

By using ^ [T ( ff o r o s \ = a s + a r we get d T q rs I TT = q r I rs + <f I Tr , so that we have 

q^ TS = \d T q rs \ rT -fft. (B12) 

By using d A d B {T*f a c a D ) = 2 Tft D we get 

2q\ rs = d 2 T q rs \ TT . (B13) 



Therefore the relevant multipoles are expressible in terms of M(i), j^, d T q rs ^ TT and 
d%.q rs \ TT , where g rs l rT " is the mass quadrupole. 

However in Eq.(6.11) we have T^ T ^ rs , whose stress tensor monopole is q( TT ^ rs = 
f d 3 aT^ T)rs (r,a), and not Tft. From Eq.(6.9) we have T^ T)rs = £ u „ T rsuv Tft. As a 
consequence we get q™ rs = J d 3 aT^ T)rs (r,a) = j d 3 a V rsuv Tft(r,a). 

By using Eq.(6.9) we get (the surface terms vanish with the assumed support and bound- 
ary conditions) 

q\" = J d 3 aTft(r,a) = J d 3 a [± 5 rs + T$ T)rs ] (r, a) = 



= q {TT)\rs + \ 5 rs ^ q \uu ^ 

3 



1 11 

-v. q(TT)\ rs — g\ rs $rs ^ q\ uu = - Q 2 (g rs I TT § rs ^ q uu I rr ) _ 

3 2 3 

u u 

(B14) 

Therefore the monopole of the TT stress tensor T^ T ^ rs is connected to the second time 
derivative of the mass quadrupole as usual. 

We also need higher multipoles q ri - rn \ rs m J d 3 a a Tl ...a r " Tft(r, a) and their TT ana- 
logues. For the dipole we get 

J d 3 ad A d B (a ri a r a s T^f{r, <r)) = <9 2 q rirslrT = 2 (V l|rs + q r ^ s + <f |rir ) . (B15) 
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We see that the second time derivative of the mass octupole is connected to a combination 
of stress tensor dipoles. For r = s = a we have Q^q r ^ aa \ TT — 2 (g r il aa 4- 2g a ' na ). However 
qn\aa canno t be expressed only in terms of the mass octupole. 

The same pattern holds for the higher multipoles q ri -- rn \ rs . In particular with same meth- 
ods used in Ref. [21] we get the following results for the stress tensor monopole g' rs and dipole 

qu\rs 

q \rs = g \sr = g r\rs^ 
q r\us + q s\ru = ^ q rs\ru ^ 

d T q rlus + d T q slur = d 2 q rslTU , 

2 (df rs + d T <f |sM + d T q slur ) = d 3 T q rsulTT , 

d T q u \ rs = -d 3 T q rsulTT + - (d 2 q urlTS + d 2 q uslrr -2d 2 <f s|ru ). 

(Bi6) 

in terms of the momentum dipole g r '™ and quadrupole q rs \ TU and mass octupole q rsu \ TT . 
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Appendix C: The Balance of Momentum and Angular Momentum for GW's 



Here we give the evaluation of the balance equations for the 3-momentum and the angular 
momentum by using the conservation of the corresponding ADM generators (4.22) and 
(4.23). 



1. The Balance of Momentum for the GW's 

By using Eqs.(6.4) and (6.8) the first term in p\ 2 )adm 111 Eq.(4.22) can be written in the 
form 

Jd 3 a ( d r R-a M- a i d T i%) (r, B) = j 'd 3 a ( £ 9 r lf> M ab lf>) (r, a) = 

ab ab 

= Jd 3 <7(Y, d r%ysdr%l s )(r,a). (CI) 



Then, like in Eq.(7.23) for the ADM energy, the ADM momentum (4.22) can be written 
in the form 

pr _ f d3(j( (matter) r { ra d) r , ^Jr \, . 

p (matter)r^ ^ = M^?(t, 5) = ^ ~ ^( r )) Vi 

i 

3 

uv 

ffii,™ = K,r , (E| r i i, -2r(')) + iE<,r )s ^.| r i' ) - 

a as 
ka(W) 1 \a{UV) \ /l (UV)dr 3 \, 

■ M l)r A M (D M (i) a K W)( T > a )- 

(C2) 



c 3 



Like after Eq.(7.23) let us divide the 3-space S T in two regions by means of a sphere S 
of big radius R » l c : a) an inner region V( inner ) with a compact sub-region V c of linear 
dimension l c containing all the particles (and the electro-magnetic field if we would add it); 
b) an asymptotic far region V(f ar y Let rf = <J r be a unit 3- vector. 

Since we have P r ADM = P^dm^ + ^adm^ and P p ™i+2) )r ( T i °)\<?ev far = 

(int) v / — *\ I n , A r V(far) r ts (rad) r i -»\ 

Pp(i+2)( T > (7 )lffe^ = 0, we get P AD { M = J v ^ ^^(r.ff). 

d P r 

Since P T ADM ~ is a constant, we have — J^ 51 = so that we get the following result in 
place of Eq. (7.24) 



90 



" ADM _ a ADM 

dr dr 

c 3 



32irG 



- d 3 ad T p r p l™ d ^{T,a) = 

d 3 a J2 [dr %l s 81 ^ rs + d r d T %l s d T %l s ] (r, a). 



(C3) 

By using the implication d T A h^ )rs = -n c d c 4 h^ rs + 0(1/ R 2 ) (n c = a c /\a\) of Eqs.(7.19), 
Eqs (7.25) are replaced by the following expression 



7 p r ^(inner) j f} T ^(f ar ) 



ADM _ J ADM 

.3 



C 



32nG 



c 3 



/ E E 9 * ( n ° d r %ls 9r % T )rs ) (r, a) + 0(1/ R 2 ) 

JV(far) rs c 



32?r ^ # / <*(«» 61) d<p J2 { d r 9r %l s ) (t, <?) + 0(l/i? 2 ) = 

7 i " ! 6 f d(cos6)d<pJ2 A s WU v(n)d T q lsw (T-R)d r q luv (T-R) + 0(l/r 2 ) 



2ttc 3 
G 

8^ 



/" d(cos9)d<p A -«» [9 T 3 ^ |rT (r-i?)] [,9 r ^g™l- (r _ jR)]+0(1/r 2 )) 

^ suini) 

(C4) 



after having used g' rs = \d 2 q rs \ TT , see Eq.(B13). This is the standard result for the momen- 
tum balance of GW's (see for instance Eqs. (1.164) and (3.83) of Ref.[21]). 

2. The Balance of Angular Momentum for GW's 

By using Eqs. (6.4) and (6.8) the last two lines of jJ^ADM m Eq.(4.23) can be written in 
the form (see Eq.(2.51) of Ref.[21], where the angular momentum is defined with an overall 
minus signa with respect to us) 

JL- j dV ( £ Br ^ (a' a, - o- a r) %n, v - 

uv 

- E (%>u d r %lu - %l u d T %^ ru )) (r, a). (C5) 

u 

Then, like in Eq.(7.23) for the ADM energy, the ADM angular momentum (4.23) can be 
written in the form 
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jrs 

J {l+2)ADM 



?3 / (matter) rs , (rad) rs , (int)rs\ i -»\ 

d a (p} (1+2) + P ;. (1+ j 2) + p) {l l 2) j (r, a), 



)](r,<?), 



E/ 4 ,TT o 4.TT _4.TT o 4> 
^ " (l)ru °V n (l)su ~ n (l)su °r It- 



tint) rs / 
Pj-(l+2) M 



A 

_r w(tfV) /p(l) _ ^ 1 M (UV) 

c 



A 



ft., dl 



(i)- A ; A c 

M C 

j E C * - 1 1 «'» + ^ | r?> - AC I r?') (r. 9). 

(C6) 



Like after Eq.(7.23) let us divide the 3-space S T in two regions by means of a sphere S 
of big radius R » l c : a) an inner region V( inner ) with a compact sub-region of linear 
dimension l c containing all the particles (and the electro-magnetic field if we would add it); 
b) an asymptotic far region V^ ar y Let n r = <7 r /\<j\ be a unit 3- vector. 



Since we have Jadm 

(int) rs 



f s V (far) 
ADM 



zrsV, 



J ADM allCL P 



(int)rs/ ->\\ n , T rs V(/ar) f ?3 (rad)rs/ ->\ 

Pj{i+2) ( r ^)\ffev far = 0, we get J ADM = J d A o p] (1+ ; 2) (r,a). 



(rad) rs , 



(matter) rs 
j(l+2) ' 



Since J^dm * s a constant, we have — 



so that we get the following result in place 



of Eq.(7.24) (again by using Eqs.(7.19) and n c = cr c /\a\) 



7 j rs V(inner) 

a J ADM 

dr 



J j rsV (far) 

u "ADM 

dr 



d ad T p j{ ± +2) >{T,(T) 

- (far) 



d 3 a 



(far) 



" C 9cP*tt (r,S) + 0(UH') 



= / 

J V (far) 



d 3 ad r 



c rs (rad) 
71 Pj{l+2) 



(r,a)+0(l/R 2 ). 



(C7) 
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As a consequence, we get 



7 j rs V(inner) 

a ADM 

dr 



dJ 



rs V (far) 

ADM 



dr 



12nG 



R 2 j d(cos9)dip[j2 (d T %^ uv (a r d s -a s d r )% 

r,a) + 0(l/R 2 ), 



TT 

)uv 



/ , y n {l)ru°T rl (l)su°T n (l)ru 



accord with Eq.(2.61) of Ref.[21] for the balance of angular momentum. 
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